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FOREWORD 


The  dense  plasma  focus  has  been  investigated  at  many  laboratories  as  a 
possible  fusion  device.  Typical  plasma  parameters  for  this  device  are  electron 
temperatures  of  1  keV,  densities  of  10*®  per  cc ,  and  confinement  times  of 
100  ns.  Characteristic  of  the  plasma  focus  discharge  are  intense  soft  and  hard 
x-ray,  optical  and  microwave  radiations.  Also  emitted  from  the  focus  are 
electrons,  ions,  and  with  certain  filling  gases,  neutrons.  The  emphasis  of  this 
work  is  to  investigate  the  electron  and  ion  emission  from  the  plasma  focus  and 
the  development  of  appropriate  diagnostics  to  accomplish  this  task.  The  view¬ 
point  of  the  focus  as  a  bipolar  diode  may  provide  a  natural  description  of  the 
observed  phenomena.  A  Mather  geometry  plasma  focus  device  has  been  constructed. 
Diagnostics  used  for  electron  measurements  include  current  monitors,  a  calori¬ 
meter,  a  magnetic  energy  analyzer,  and  an  emittance  meter.  Ion  beam  measure¬ 
ments  were  performed  using  a  fluxmeter,  a  pinhole  camera,  an  ion  beam  emittance 
meter,  and  a  Thomson  parabola  analyzer.  In  addition,  opening  switch  behavior  has 
been  investigated.  Intense  electron  beams  are  observed  with  pulse  lengths 
approximately  5  ns  at  energies  of  20  to  200  keV  and  currents  of  1  kA  and  above. 
The  root-mean-square  (rms)  emittance  of  the  electron  beam  is  found  to  be 
approximately  1  cm-rad.  Ion  measurements  show  copious  ion  emission  of  very 
energetic  ions  of  up  to  several  MeV.  A  measurement  of  rms  emittance  of  the  ion 
beam  was  made.  Results  support  the  bi-polar  diode  model  of  the  plasma  focus 
device . 


I  would  especially  like  to  thank  Dr.  Moon-Jhong  Rhee  for  his  guidance  and 
direction  throughout  this  experimental  work.  Also,  I  would  like  to  thank 
Dr.  John  R.  Smith,  Dr.  Won  Namkung,  and  Mr.  Cheng-Mu  Luo  for  their  collaboration 
while  working  on  these  experiments. 
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CHAPTER  1 

OVERVIEW  OF  PLASMA  FOCUS  RESEARCH 


INTRODUCTION 

The  plasma  focus,  developed  in  the  Soviet  Union  and  the  United  States  by 

1  o 

Filippov  and  Mather,  is  basically  a  linear  pinch  in  a  geometry  which  allows 

rapid  compression  of  the  plasma  prior  to  the  pinch  by  magnetic  self  forces.  In 

addition  to  the  high  plasma  densities  of  10^^  per  cc^’^'^  and  temperatures  of  a 

few  keV  ’  which  may  last  a  hundred  ns,  many  interesting  phenomena  are 

789  8  10 

observed.  Among  these  are  intense  radiations  in  the  microwave,  ’  ’  optical,  ' 
soft,®’^^  and  hard^’^’®’^^'^^  x-ray  regions.  Electron,^ ’ ^ 

ion,3"6,16,18-21,25-34  ^^^^^^,^^.^3-7 , 16 , 18 , 21 , 30-32 , 35  emissions  have  also  been 


investigated  at  many  laboratories. 


The  plasma  focus  originated  from  a  simple  coaxial  plasma  accelerator 

concept.  Mather  investigated  a  higher  pressure  mode  of  operation  of  the 

device^  while  Filippov's  investigation  centered  around  new  methods  of  breakdown 

across  insulators  for  z-pinch  applications.^  When  it  was  observed  that  copious 

neutrons  were  produced  by  this  device,  investigations  began  into  possible  use  of 

this  device  for  fusion  applications;  however,  it  was  soon  realized  that  an  ion 

beam- target  interaction  explained  the  neutron  distribution  better  than  a 

3  18  31 

thermonuclear  model.  ’’ 


The  plasma  focus  discharge  consists  of  three  main  stages:  (!)•  breakdown 

across  the  insulator,  (2)  rundown  stage,  and  (3)  collapse  to  focus  or 

instability.  The  breakdown  is  initiated  across  the  insulator  in  Figure  1-1 

through  a  low  pressure  gas  (0.1 -5.0  Torr)  by  closing  a  switch  connected  to  a  high 

voltage  capacitor.  It  has  been  observed  that  in  some  modes  of  operation  the 

breakdown  may  occur  between  the  inner  and  outer  electrodes  away  from  the 

37  39 

insulator  surface  and  may  show  some  filamentary  structure.  This  regime. 
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which  may  depend  on  filling  pressure  and  insulator  geometry,  is  associated  with 
poor  focus  formation  and  should  be  avoided.  It  has  been  suggested  by  Krompholz, 
et  al . ,  that  a  circular  knife  edge  placed  at  the  junction  of  the  outer 
electrode  and  the  base  of  the  insulator  may  serve  to  enhance  the  E- field  and  lead 
to  uniform  breakdown  which  will  improve  the  final  focus  formation 
reproducibility. 


The  rundown  stage  occurs  as  the  current  sheath  lifts  of  the  insulator  surface 

and  is  driven  down  the  center  electrode  length  by  forces.  Tn  the  snowplow 

model,  as  the  current  sheath  advances  through  the  gun  region,  all  the  gas  in  the 

gun  region  is  pushed  out  leaving  vacuum  behind.  The  current  sheath  travels  at  an 

approximately  constant  velocity  of  a  few  cm//isec  which  is  determined  by  momentum 

balance  between  the  filling  gas  and  the  magnetic  pressure  behind  the  current 

sheath.  In  cylindrical  coordinates,  the  parabolic  shape  of  the  current  sheath  is 

39 

given  by  the  function  : 


z(r)-  - 


2i 


k(k2  .  i2)V2.  1,  kLUli;  , 

l-t-(l-i^)/2 


(1-1) 


where  a  is  the  center  electrode  radius,  k  is  r/a,  and  i  is  an  experimentally 
determined  parameter  given  by  the  momentum  balance  condition: 


i 


o  2  2  2 

8njr  a  v 


(1-2) 


where  I  is  the  current  flowing  through  the  sheath,  n  is  the  density  of  gas  and  v 
is  the  velocity  of  the  sheath.  The  quantity  i  approximately  equals  unity. 

The  m-0  instability,  which  has  the  highest  growth  rate,^*^  occurs  after  the 

current  column  carrying  100  kA  to  1  MA  collapses  toward  the  axis  due  to  J  B. 

z  6 

forces.  It  is  during  this  instability  stage  that  most  physically  interesting 

19  3 

phenomena  occur.  From  the  hot  (T^-l  keV) ,  dense  (n-10  cm  )  plasma,  intense 
optical,  soft  x-ray  and  microwave  radiations  are  observed.  In  addition, 
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electrons  and  ions  are  accelerated  in  the  focus  to  energies  of  a  few  hundred 
keV  per  charge  which  have  been  attributed  to  the  m-0  instability  by  some 
models. Byproducts  of  these  are  hard  x-rays  from  brehmstrahlung  as  the 
electrons  strike  the  copper  electrodes;  also,  (if  the  filling  gas  is  deuterium) 

O  Q 

10-10  neutrons  per  shot  are  produced  by  ion  beam-target  spallation  reactions. 
Semi -empirical  models  of  the  m-0  instability  in  the  plasma  focus  by  R.  Deutsch^^ 
and  Hohl  and  Gary^^  indicate  that  a  voltage  on  the  order  of  100  kV  may  be  induced 
in  the  focus  by  the  m-0  instability.  These  models  assume  some  typical  m-0 
constriction  velocities  observed  with  streak  camera  photography  or  current 
reduction. 

It  has  been  observed  with  laser  scattering  experiments^^  that  subsequent  to 
the  onset  of  the  m-0  instability,  there  is  a  rapid  increase  in  the  electronic 
temperature  accompanied  by  a  period  of  resistive  magnetic  field  diffusion  through 
the  plasma  to  the  axis.  It  has  been  speculated  that  lower  hybrid,  or  electron 
cyclotron  drift  instabilities  may  be  responsible  for  some  of  the  heating  '  and 
subsequently,  ion  acoustic  instability  may  dramatically  increase  the  resistivity 
of  the  plasma.  In  fact,  recently  microwave  emission  has  been  observed  in  the 
plasma  focus  which  may  be  attributed  to  the  lower  hybrid  drift  instability.^® 
Figure  1-2  shows  a  temporal  progression  of  these  phenomena. 

Particle  energies  of  up  to  several  MeV  have  been  observed  but  the  yields 
at  these  energies  are  typically  a  few  orders  of  magnitude  below  those  in  the 
100  keV  range.  The  physical  mechanism  of  the  intense  particle  beam  generation  in 
a  plasma  focus  device  is  unclear  up  to  the  present  time.  The  goal  of  this  work 
is  to  investigate  the  physical  nature  of  the  plasma  focus  electron  and  ion  beams 
in  hopes  that  an  understanding  of  this  diode-like  behavior  will  lead  to  an 
improved  physical  understanding  of  the  plasma  focus  device  as  an  excellent  source 
of  electron  and  ion  beams  in  the  regime  considered.  This  requires  that  the 
source  of  the  most  intense  beams  be  investigated. 

PLASMA  FOCUS  ACCELERATOR 

A  plasma  focus  device  may  be  described  by  a  simple  lumped  circuit  model  as  in 
Eltgroth^^^  and  Mather.^®  In  Figure  1-3,  Cgj.^  is  the  external  capacitance  which 
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consists  mostly  of  the  capacitance  of  the  main  capacitor  bank.  Stray  capacitance 
is  negligibly  small.  is  the  external  circuit  inductance  which  arises  from 

the  bank  inductance  and  all  other  circuit  inductances  exclusive  of  the  gun 
region-plasma  sheath  inductance.  Rgxt  external  circuit  resistance 

excluding  the  gun  region-plasma  sheath  resistance.  Lp  is  the  gun-plasma  sheath 
inductance  which  is  a  variable  quantity,  changing  with  the  geometry  of  the 
current  sheath.  Rp  is  the  gun-plasma  sheath  resistance.  Finally,  R^  is  either 
the  diode  load  resistance  or  the  leakage  resistance  across  the  insulator  which  is 
in  parallel  with  Rp.  At  time  t-0  the  switch  is  closed  and  the  current  begins  to 
flow  through  the  circuit.  Subsequent  to  the  current  sheath  breakdown  across  the 
insulator  in  the  gun  region,  Rp  is  on  the  order  of  milliohms  (Spitzer 
resistivity)  as  pointed  out  by  Mather.  Rgxt’  spark  gap 

switch,  might  be  much  higher.  Lp  is  calculable  and  found  to  be  less  than  or  the 
same  order  of  magnitude  as  ,  depending  on  the  position  and  shape  of  the 

current  sheath.  is  a  dominant  element  by  the  nature  of  the  device.  The  net 

effect  of  these  circuit  elements  is  to  produce  a  ringing  circuit  which  is 
underdamped  by  the  total  reistance  which  may  mostly  appear  in  the  spark  gap  and 
the  plasma  current  sheath.  The  instability  occurs  as  the  sheath  collapses  toward 
the  axis  at  which  time  the  plasma  resistivity  rapidly  rises  and  the  current 
flowing  through  the  circuit  begins  to  drop.  A  significant  induced  voltage 
appears  across  the  diode  load  Rj^  at  this  time.  It  may  be  shown  that  in  the  focus 
stage  L  ^  is  much  larger  than  I  In  fact,  it  is  important  to  note  that  the 

signs  of  these  two  terms  are  opposite  each  other  during  the  plasma  pinch 
stage.  I  ^  opposes  the  flow  of  current^^  and  L  ^  favors  production  of  the 
beam.  Another  important  point  to  be  considered  is  the  leakage  current  across  the 
insulator.  The  induced  voltage  that  appears  at  the  time  of  the  focus  will  also 
appear  across  the  insulator.  If  the  current  sheath  has  not  swept  out  all  of  the 
gas  between  the  electrodes,  there  may  be  a  possibility  of  a  subsequent  breakdown 
across  the  insulator  which  will  immediately  short  circuit  the  diode  load. 

If  the  m-0  instability  completely  interrupts  the  current  flow  at  the  time  of 
the  instability  the  voltage  approaches  infinitely  large  in  this  model.  The 
circuit  model  requires  some  further  input  to  describe  the  form  of  the  plasma 
resistance  as  a  function  of  time  in  order  to  predict  the  diode  voltage  and 
current.  Unfortunately,  the  theory  of  the  plasma  focus  instability  is  not 
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sufficient  at  the  present  time  to  accomplish  this  task.  This  is  the  motivation 
for  the  development  of  a  transmission  line  theory  of  the  plasma  focus. 

A  transmission  line  theory  of  the  plasma  focus  was  developed  in  an  attempt  to 
explain  the  presence  of  a  finite  voltage  pulse. In  essence,  a  transmissionline 
theory  accounts  for  the  presence  of  an  electromagnetic  field  and  the  dissipation 
of  its  energy  in  a  resistive  load.  The  simplest  possible  model  for  the  plasma 
focus  is  found  in  Figure  1-4.  Immediately  after  the  formation  of  the  current 
sheath,  the  plasma  focus  is  modeled  as  a  current  charged  transmission  line  shown 
in  Figure  1-4  with  the  switch  at  1  closed.  This  assumes  that  the  resistance  of 
the  current  sheath  and  sparkgap  are  negligible,  and  all  characteristic  impedances 
of  each  section  of  the  plasma  focus  device  are  matched  to  Zq .  At  this  stage,  the 
current  charged  transmission  line  can  be  thought  of  as  two  oppositely  traveling 
waves  of  current: 

I  -  I_^(t-  J)  +  I_(t+  ,  (1-3) 

where  the  magnitude  of  and  I_  is  Iq/2  (half  the  total  charging  current),  and  v 
is  the  velocity  of  transverse  electromagnetic  (TEM)  wave  propagation  in  the 
transmission  line.  The  voltage  waves  are  of  the  form: 


V  -  I 


-  '-/"c 


(1-4) 


where  Zg  -  7  (1^/C^) ,  1^  is  the  inductance  per  unit  length,  and  is  the 
capacitance  per  unit  length.  At  time  t-0,  the  switch  is  opened.  The  reflection 
and  transmission  coefficients  determined  by  boundary  conditions  in  this  simple 
model  are: 


and 


(1-5) 


%  -  (R^  +  Zq)  ’ 

(satisfying  the  boundary  conditions  1  +  and  conservation  of  energy 

2  2 

“  1):  hence,  if  the  load  resistance  Rj^  is  matched  to  the  impedance  of  the 
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transmission  line,  equals  0  and  equals  1.  Subsequently,  the  right 
traveling  current  wave  in  Figure  1-4  is  dissipated  in  the  load  resistance.  The 
output  pulse  is  of  duration  21/v,  voltage  IqZq/2,  and  current  1q/2.  It  is 
important  to  emphasize  that  these  parameters  arise  from  an  ideal  opening  switch 
and  perfect  impedance  matching.  These  depend  upon  the  detailed  physics  of  the 
focus  instability  which  are  not  well  understood. 

As  a  simple  non- ideal  case,  we  may  suppose  that  the  load  resistance  is  not 
precisely  matched  to  that  of  the  transmission  line.  In  this  case,  the  reflection 
and  transmission  coefficients  are  no  longer  0  and  1,  respectively.  The  two 
distinctly  different  regimes  are  >  Zq  and  Rl  <  Zq.  If  Rj^  is  greater  than  Zq, 
then  p^  is  greater  than  0  and  is  greater  than  1.  It  can  easily  be  seen  that 
the  output  voltage  will  consist  of  multiple  pulses  of  changing  polarity,  the  i^^ 
pulse  having  amplitude  IqZq  (-p^)^  ^  i^+P^)/2  and  duration  21/v  as  depicted  in 
Figure  1-5.  On  the  other  hand,  if  R^^  <  Zq,  p  is  less  than  0  and  the  output  will 
all  be  of  one  polarity.  See  Figure  1-5.  The  distribution  of  this  current  output 
between  electron  and  ion  currents  will  depend  upon  the  diode  characteristics.  We 
will  describe  the  bi -polar  diode  in  the  following  section. 


The  acceleration  of  ions  and  electrons  across  a  diode  has  been  investigated 

quite  extensively.  The  simple  relativistic  diode  model  for  the  plasma  focus 

maybe  a  one - dimens ional ,  bi -polar,  Child- Langmuir  diode.  In  this  model,  two 

semi-  infinite  slabs  of  plasma  are  separated  by  a  fixed  distance  with  a  fixed 

voltge  supplied  between  them.  The  resultant  current  flow  is  space  charge 

52 

limited.  The  ratio  of  ion  current  to  electron  current  is  given  by; 


^ion  ,  ^  (7  1)  1V2 

I  "  ^  M.  2  ‘ 

e  1 


(1-7) 


where  m^  and  are  the  electron  and  ion  rest  mass,  (7  +  1)  m^c  /e  is  the  voltage 
across  the  diode,  c  is  the  velocity  of  light  in  vacuum,  and  e  is  the  electronic 
charge . 


A  correction  including  two-dimensional  effects  has  been  proposed  by  Goldstein 
and  Lee^  and  applied  to  the  plasma  focus  by  Gullickson.  In  the  so-called 
pinch-reflex  diode,  pinching  of  the  electron  beam  across  the  gap  enhances  ion 
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current  somewhat,  but  electron  current  still  dominates  the  net  current  flow 


across  the  diode  for  most  geometries 

2  2  V? 

current  is  the  same  as  above  multiplied  by  the  factor  (r^  +  ^ 


The  ratio  of  ion  current  to  electron 

where  r. 


IS 


the  radius  of  the  cylindrical  diode  (plasma)  and  d^j  is  the  gap  spacing.  This 
factor  arises  from  the  geometry.  Due  to  the  pinching  of  the  electron  flow,  it 
has  an  effectively  longer  gap  to  traverse  than  the  ions.  Unfortunately,  it  is 
quite  difficult  to  determine  the  interrelationships  between  r^,  d^,  7,  and 
currents  in  the  real  plasma  focus  diode  in  which  these  quantities  are  not 
fixed.  It  is  hoped  that  these  experimentally  observed  quantities  will  help  in 
the  eventual  physical  understanding  of  the  plasma  focus  diode  in  order  that  the 
very  attractive  features  found  in  this  compact  accelerator  may  find  many 
applications . 


THE  EXPERIMENTAL  DEVICE 


A  prototype  plasma  focus  accelerator  has  been  constructed  in  order  to 
investigate  the  diode -like  behavior  of  the  instability.  The  major  systems  of 
this  accelerator  are:  (1)  a  capacitor  bank  and  charging  system,  (2)  a  triggering 
system,  (3)  a  transmission  line  section,  (4)  a  diagnostics  section,  and  (5)  a 
Mather  geometry  plasma  gun.  See  Figure  1-6. 

The  capacitor  bank  used  is  either  a  single  3  kJ ,  20  kV  low  inductance 
(40  nH)  capacitor  or  a  set  of  four  connected  in  parallel.  This  bank  is  charged 
through  a  500  Kilo  Ohm  resistor  by  a  Universal  Voltronics  model  bal  22-35  power 
supply.  The  device  is  triggered  by  a  spark  gap  switch  with  a  midplane  trigger 
electrode  shown  in  cross  section  in  Figure  1-7.  Typically  the  gap  utilizes 
nitrogen  gas  at  or  slightly  above  atmospheric  pressure.  The  breakdown  field  is 
given  by: 

E^^^(kV/cm)  -  [24.6p  +  6 . 7  (p/d^^^)^/2  1  (1.8) 

where  p  is  the  pressure  in  atmospheres,  is  the  effective  electrode  gap 

separation  and  F  is  the  field  enhancement  factor,  i.e.,  the  ratio  of  the  maximum 
to  the  mean  field  on  the  electrodes.  For  our  spherical  electrodes,  is 

0.115d,  where  d  is  the  actual  gap  spacing,  and  the  field  enhancement  factor  F  is 
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1.8.^^  The  gap  spacing  used  varied  somewhat,  however  remained  near  1  cm  for  most 
experiments.  An  estimate  of  switch  inductance  given  by  Miller  gives  14  nH/cm 
hence  our  switch  adds  about  14  nH  to  the  total  inductance.  A  Thyratron  tube 
circuit  is  used  to  provide  the  high  voltage  (18-24  kV)  pulse  to  the  trigger 
electrode  to  initiate  gap  breakdown. 

The  transmission  line  section  consists  of  a  30  cm  long  coaxial  line  with  an 
inner  conductor  diameter  of  5.1  cm  and  an  outer  conductor  diameter  of  14.6  cm. 

The  total  inductance  is  64  nH  for  this  section.  The  characteristic  impedance 
(Zq)  is  ^  ohms  for  this  geometry,  where  e  is  the  dielectric  constant  of  the 
medixim  between  the  conductors.  For  a  water  dielectric,  this  gives  Zq  -  7  ohms, 
while  for  air,  60  ohms.  Water  is  used  exclusively  in  this  work  due  to  is 
conductive  properties  assisting  in  holding  off  a  corona  discharge  across  the 
insulator  before  firing,  therefore  the  matched  output  pulse  would  be  18  nsec  in 
duration  and  the  output  voltage  would  be  one  half  the  charging  current  times  7 
ohms . 

A  2.5  cm  long  diagnostics  section  is  inserted  between  the  transmission  line 

56-59 

and  the  coaxial  gun.  See  Figure  1-8.  A  Rogowski  coil  is  made  from 

2 

30  windings  on  a  ring  with  a  9  mm  cross  section  and  a  mean  radius  of  7.65  cm. 

The  coil  signal  is  integrated  with  a  50  ^is  (-  measured  R,C)  passive  integrating 
circuit.  The  calibration  is  calculated  to  be  68.55  kA/V.  A  capacitive  voltage 
divider  probe  (D-dot)  is  constructed  by  placing  a  1  cm  radius  disc  1  cm  from  the 
inner  surface  of  the  diagnostics  ring.  The  signal  from  the  probe  is  integrated 
by  a  passive  2.63  /is  integrator.  The  calibration  is  calculated  to  be  19.2  kV/V 
with  water  dielectric  in  the  transmission  line.  The  B-dot  probe  was  not  used  in 
any  experiments  discussed  in  this  work. 

The  usual  Mather  geometry  plasma  gun  is  shown  in  Figure  1-9.  It  consists  of 
a  7.5  cm  inner  diameter  copper  outer  electrode  (cathode)  and  a  2. 5  cm  diameter 
inner  electrode  (anode).  For  most  experiments,  a  squirrel  cage  outer  electrode 
is  used  which  is  made  from  12,  0.125 -inch  copper  rods  spaced  equally  on  a  7.5  cm 
diameter.  A  pyrex  glass  insulator  separates  the  inner  electrode  from  the  brass 
backplate  which  is  electrically  the  cathode.  The  insulator  is  fixed  in  place  on 
the  plexiglass  with  silicone  sealant.  The  2.6  cm  long  insulator  surface  is  made 
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flush  with  the  cylindrical  surface  of  the  inner  electrode  with  a  small  gap 
between  them  (1  mm)  to  prevent  mechanical  shock  on  the  inner  electrode  from 
shattering  the  glass.  This  geometry  is  utilized  for  ion  beam  measurements  since 
the  ions  are  accelerated  away  from  the  anode.  In  order  to  investigate  the 
electron  beam  which  is  accelerated  down  the  middle  of  the  hollow  center 
electrode,  we  must  invert  the  usual  gun  geometry  as  shown  in  Figure  1-10.  In 
this  experimental  arrangement,  the  capacitor  bank  is  charged  negatively  so  that 
the  center  conductor  on  the  transmission  line  becomes  the  plasma  focus  gun 
cathode . 
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FIGURE  1- 


1)  BREAKDOWN 


2)  RUNDOWN 


3)  FOCUS 


MATHER  GEOMETRY  PLASMA  FOCUS  GUN  SHOWING  BREAKDOWN, 
RUNDOWN  AND  FOCUS  STAGES  OF  CURRENT  SHEATH 
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1)  CURRENT  COLUMN 


ANODE 


2)m  =  0,  ACCELERATION 


3)  MICROINSTABILITY,  ACC. 


FIGURE  1-2.  SEQUENCE  OF  INSTABILITY  IN  THE  PLASMA  FOCUS 
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FIGURE  1-3.  LUMPED  CIRCUIT  MODEL  OF  THE  PLASMA  FOCUS.  AS  IHc  t-LASMA 
RESISTANCE  INCREASES,  A  HIGH  VOLTAGE  APPEARS  ACROSS  THE 
DIODE  LOAD  AND  THE  INSULATOR. 
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FIGURE  1-6.  PLASMA  FOCUS  DEVICE 
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ROGOWSKI  COIL 


D-DOT  PROBE 


FIGURE  1-9.  DETAIL  OF  THE  MATHER  GEOMETRY  GUN 
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CHAPTER  2 

OPERATION  CHARACTERISTICS  OF  THE  PLASMA  FOCUS 


BREAKDOWN  STAGE 

A  symmetric  and  uniform  initial  breakdown  across  the  pyrex  insulator  in  the 
plasma  gun  may  be  quite  important  to  the  formation  of  a  symmetric  and 
reproducible  focus.  The  breakdown  occurs  as  a  high  voltage  appears  between  the 
anode  and  cathode  subsequent  to  the  triggering  of  the  spark  gap  switch.  The 
conditioning  process  of  the  insulator  has  been  evident  in  many  experiments. 
Typically,  a  new  plasma  focus  device  insulator  requires  anywhere  from  0  to  60 
shots  before  a  reproducible  current  drop  in  the  main  transmission  line  Rogowski 
coil,  or  a  corresponding  voltage  spike,  is  evident.  This  process  of  condition¬ 
ing  has  been  observed  in  vacuum  spark  gaps;^®  however,  we  wish  to  include  the 
plasma  focus  device  as  one  undergoing  the  similar  type  of  conditioning  process. 

In  the  plasma  focus  device,  high  currents  well  above  100  kA  are  passed  between 
the  electrodes  through  a  very  low  pressure  (sub  Torr)  gas.  During  this  time, 
heating  and  sputtering  of  the  electrode  surfaces  would  result  in  the  formation  of 
a  metallic  vapor  in  the  gun  region  which  would  preferentially  redeposit  on  the 
electrodes  due  to  their  cooler  temperatures  (since  their  thermal  conductivity  is 
higher  than  that  of  the  pyrex  insulator).  However,  eventually  some  electrode 
material  would  be  deposited  on  the  pyrex  insulator,  and  subsequent  shots  would 
then  preferentially  deposit  on  the  regions  which  already  have  some  metallic 
deposit.  This  process  would  lead  to  a  clumping  of  the  metallic  deposits  (or 
perhaps  crystal  growth)  in  these  regions. 

It  was  observed  that  our  pyrex  insulator  became  coated  with  a  non-conducting 
metallic  copper  layer  in  our  electron  beam  geometry  gun  which  has  all  copper 
electrodes.  This  gun  worked  very  reproducibly  for  several  hundred  shots  before 
the  insulator  broke  due  to  mechanical  shock.  When  a  chip  of  this  insulator  was 
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placed  under  an  optical  microscope,  it  was  noticed  that  the  surface  showed 
crazing,  in  addition  to  the  coating  of  copper.  A  picture  of  the  insulator  taken 
through  an  optical  microscope  is  shown  in  Figure  2-1. 

2 

Regions  of  approximately  0. 1x0.1  mm  were  evidently  separated  by  the  small 
surface  cracks  which  may  have  served  to  maintain  the  overall  non- conductive 
property  of  the  insulatorsurface .  Under  the  scanning  electron  microscope, 
further  detail  of  the  coating  is  evident  in  Figure  2-2,  indicating  that  there  are 
many  micron-size  clumps  of  metal  which  seem  to  be  preferentially  deposited  along 
the  surface  cracks . 

It  is  hypothesized  that  these  regions  of  metallic  coating  may  serve  to  grade 
the  insulator  naturally,  providing  a  uniform  and  symmetric  breakdown  along  the 
surface  which  leads  to  a  symmetric  and  reproducible  focus  formation. 


CURRENT  SHEATH  RUNDOWN  STAGE 


The  sheath  rundown  stage  is  investigated  with  the  use  of  a  framing  camera 
which  has  an  exposure  time  of  10  nsec  and  a  frame  separation  of  50  nsec.  As 
pointed  out  previously,  the  velocity  of  the  current  sheath  is  mainly  determined 
by  momentum  balance : 


V  -  ( 


2  2 
8njr  a 


ii 


4.2  cm//js , 


(2-1) 


-  3  -  3 

where  the  values  I  -  205  kA,  n  -  2.35  x  10  kg/m  ,  and  a  -  1.27  cm  have  been 
used  for  the  1  Torr  argon  (Ar)  shot  which  we  shall  analyze.  The  average  velocity 
of  the  sheath  is  determined  during  three  different  regimes  in  the  shot  from 
photographs  similar  to  those  shown  in  Figure  2-4.  The  three  regimes  are:  (1) 
radial  expansion  outward  from  the  insulator  (1.5-2. 5  /is),  (2)  motion  down  the 

coaxial  electrodes  (2. 5-3.4  ps  ),  and  (3)  radial  collapse  at  the  tip  of  the 
center  electrode  (2.9-3. 5  /is).  In  regimes  (1)  and  (3),  the  image  converter 

camera  is  placed  on  axis  of  the  system  so  that  only  the  radial  component  of  the 
velocity  is  measured  as  in  Figure  2-3(A).  For  regime  (2),  the  camera  is  placed 
so  that  a  side-on  view  of  the  sheath  is  obtained  as  in  Figure  2-3(B).  The 
average  velocity  of  the  radial  expansion  from  1.5  to  2.5  ^s  is  found  to  be 


2-2 


NSWC  TR  88-394 


1.5  (+/-  0.2)  cm//is  (the  velocity  along  the  center  conductor),  and  it  is  found 
to  average  1.4  (+/-0.2)  cm/ps  from  2.5  to  3.4  ps.  Finally,  the  radial  collapse 
is  found  to  average  2.6  (+/-  0.2)  cm/ps  during  the  time  from  2.9  to  3.5  ps . 

The  error  limits  are  approximately  determined  from  the  re -solution  of  the  framing 
photographs.  The  value  found  for  case  (2)  shows  approximate  agreement  with  the 
simple  existing  theory  which  gives  4.2  cm/ps. 

OPENING  SWITCH  STUDIES 

The  plasma  focus  device  may  be  viewed  as  an  opening  switch  connected  to  a 
load  (i.e.,  plasma  diode). Characteristics  of  the  opening  switch  may  be 
investigated  with  the  transmission  line  Rogowski  and  voltage  probes.  A  schematic 
of  the  plasma  focus  is  found  in  Figure  2-4.  The  device  is  represented  by  a 
current  charged  transmission  line  with  impedance  Zq .  Initially,  the  m-0  switch 
is  closed,  and  the  transmission  line  will  charge  to  a  current  peak  through  LC 
ringing  determined  by  the  lumped  circuit  inductance  and  capacitance  of  the 
transmission  line  and  capacitor.  When  the  current  flowing  through  the  plasma 
(closed  m-0  switch)  reaches  its  peak  (Io)t  the  m-0  instability  opens  the  switch 
forcing  the  dissipation  of  the  stored  electromagnetic  field  energy  through  the 
load  R^  (plasma  diode)  resulting  in  the  acceleration  of  ions  and  electrons  in  the 
plus  and  minus  z  directions. 

The  characteristic  impedance  of  the  transmission  line  is  found  from  the 
output  voltage  and  current  drop  (see  Figure  2-5).  The  slope  of  the  graph  from 
aleast  squares  fit  is  6.88  ohms.  For  comparison  the  calculated  value  of 
impedance  y(L^/C^)  is  7  ohms. 

The  diode  load  resistance  may  be  found  from  the  output  current  and 
voltage.  For  example,  the  shot  in  Figure  2-6  shows  an  output  voltage  of  150  kV 
with  an  output  current  of  157  kA  giving  Rl  “  0.96  ohms.  Figure  2-6  also  shows 
the  theoretical  current  and  voltage  waveforms  assuming  an  ideal  transmission  line 
behavior  with  these  values.  The  difference  indicates  a  non- ideal  behavior  which 
may  be  attributed  to  the  non- ideal  opening  switch  action  of  the  plasma  focus. 
Indeed  the  voltage  waveforms  in  Figure  2-7  indicate  a  wide  shot-to-shot 
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variation,  sometimes  showing  multiple  peaks  as  well  as  single  voltage  peaks, 
notin  agreement  with  the  simple  transmission  line  theory  except  that  the  pulse 
width  is  of  the  same  order  of  magnitude.  In  a  typical  series  of  21  shots  of 
He-Armixture ,  about  38  percent  of  the  shots  showed  single  spikes,  43  percent 
showed  double  spikes,  and  19  percent  showed  more  than  two  peaks.  The  average 
FVHM  of  a  single  spike  is  21  ns  with  a  5  ns  standard  deviation.  The  average  peak 
to  peak  separation  of  double  spikes  is  17  ns  with  a  standard  deviation  of  6  ns . 
This  erratic  behavior  is  one  of  the  most  important  contributing  factors  in  making 
the  exact  acceleration  mechanism  difficult  to  determine. 
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FIGURE  2-1.  OPTICAL  MICROSCOPE  VIEW  OF  INSULATOR  SURFACE 


FIGURE  2-2.  SCANNING  ELECTRON  MICROSCOPE  VIEW  OF  INSULATOR  SURFACE 
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lA)  TAKEN  ALONG  THE  Z-AXIS 

Ar  1.0  TORR.  DELAY  TIME  3.0  ms 


SEQUENCE 

2  4  6  8  10  12 
1  3  5  7  9  11 


(B)  SIDE-ON  VIEW 


He  1.0  TORR.  DELAY  TIME  2.8  ms 

FIGURE  2-3.  FRAMING  PHOTOGRAPHS  OF  THE  PLASMA  FOCUS 
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TIC  OF  THE  PLASMA  FOCUS 


Output  Vol tage  ( kV ) 


FIGURE  2- 


NSWC  TR  88-394 


Current  Drop  (kA) 


5.  PLOT  OF  OUTPUT  VOLTAGE  VERSUS  CURRENT  DROP  FOR  THE  PLASMA 
FOCUS  DEVICE  TO  FIND  THE  TRANSMISSION  LINE  IMPEDANCE 
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(A)  TRANSMISSION 
LINE  CURRENT 


(B)  MODEL  TRACES 


lo  =  178  kA 


1  AU/div 


TRANSMISSION 
LINE  VOLTAGE 


20  nsec/div 


Zo  =  V/AI  =  7.3  n 


RL  “  -  0.96  ^2 


FIGURE  2-6.  ANALYSIS  FOR  DIODE  LOAD  RESISTANCE  R^  (COMPARE  WITH 
FIGURE  1-5) 
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TRANSMISSION  LINE  PROBES 


990  mTorr  He 
+  10  mTorr  Ar 


500  mTorr  He 
+  500  mTorr  Ar 


68.56  kA/Div  38.4  kV/Div 

1  fisec/div  50  or  100  ns 

FIGURE  2-7.  TRANSMISSION  LINE  DIAGNOSTICS  SIGNALS  (INTEGRATED).  CURRENT 
IS  ON  THE  LEFT  AND  VOLTAGE  ON  THE  RIGHT. 
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CHAPTER  3 

ELECTRON  BEAM  MEASUREMENTS 


ELECTRON  BEAM  CURRENT  MEASUREMENTS 


A  Rogowski  coil  has  been  used  to  measure  the  current  associated  with  the 

22 

electron  beam  produced  by  the  plasma  focus  device.  The  coil  consists  of  12 
windings  on  a  plexiglas  (poly-methyl  methacrylate)  ring  of  3  mm  x  3  mm  cross 
section  with  a  3.5  cm  mean  diameter.  A  voltage  is  induced  between  the  ends  of 
the  windings  by  the  time  varying  magnetic  field  produced  by  the  current  pulse 
which  passes  through  the  ring.  The  induced  voltage  is: 


A  ^ 
r  dt 


(3-1) 


where  A  is  the  cross-sectional  area  of  the  ring,  N  is  the  number  of  windings, 

/iQ  is  the  permeability  of  free  space,  and  r  is  the  mean  radius.  This  signal  is 
passed  through  a  simple  passive  integrating  circuit.  Exact  analysis  of  the 
circuit  shows  the  output  voltage  will  be  given  by: 


V  1/RC  (f  V.  dt  -  r  V  ^dt)  .  (3-2) 

out  ■'in  out 

It  is  readily  apparent  that  if  the  current  pulse  timescale  is  much  shorter  than 
the  RC  time  constant,  the  second  term  in  the  above  equation  is  negligible  and  the 
Rogowski  output  will  be  proportional  to  the  net  current  flowing  through  the 
coil.  As  the  pulse  time  approaches  RC,  the  correction  term  is  no  longer 
negligible . 

1  ^  ^0  A 

With  an  RC-875  nsec  integrator,  the  proportionality  constant,  r— —  — r -  -  . 

RC  2  TT  r 

equals  729.1  amps/volt.  The  actual  calibration  is  found  to  be  709  amps/volt  by 
passing  a  known  current  pulse  through  the  coil. 
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The  Rogowski  coil  is  placed  directly  at  the  output  of  the  hollow  center 
electrode  drift  tube  as  shown  in  Figure  1-10. 

The  signal  risetime  is  calculated  to  be  0.3  ns,  which  is  adequate  for  the 
e-beam.  Sample  results  are  shown  in  Figure  3-1.  Figure  3-2  is  a  plot  of  the 
peak  current  as  a  function  of  gas  pressure  for  helium,  nitrogen,  and  argon 
gas.  It  should  be  noted  that  due  to  the  pressure  dependence  of  the  current 
sheath  velocity,  the  focus  will  occur  earlier,  i.e.,  at  lower  plasma  currents 
(Ip£)  at  lower  pressures.  This  data  is  plotted  in  Figure  3-3.  One  interesting 
feature  of  the  current  waveforms  is  the  long  time  scale  decay  which  indicates 
the  flowing  of  plasma  current  through  the  Rogowski  coil.  As  previously 
described,  the  current  decay  time  read  from  this  diagnostic  requires  correction 
due  to  the  finite  RC  time  of  the  integrator;  however,  the  peak  value  of  current 
which  occurs  quite  early  in  the  pulse  may  be  read  quite  accurately  without 
correction . 

A  Faraday  cup  may  provide  a  more  sensitive  current -measuring  device 

suitable  for  measuring  the  very  short  time  scale  beam  current  (10  ns).  A 

6.7  cm  diameter,  1  cm  thick  graphite  block  is  connected  to  a  0.05024  ohm  T-+M 

Research  Products  current  viewing  resistor.  In  order  to  maintain  high  vacuum 

(10*^  Torr)  in  the  Faraday  cup  region,  a  foil  is  placed  2  cm  downstream  of  the 

center  electrode  to  isolate  it  from  the  plasma  focus  static  filling  gas.  The 

Faraday  cup  is  placed  1  cm  downstream  of  the  foil.  This  serves  to  insure  that 

no  beam  filling  gas  interaction  be  recorded  at  the  Faraday  cup.  The  foils  may 

also  act  as  filters  in  that  the  foil  thickness  is  comparable  to  the  ranges  of 

electrons  in  the  foil  material.  In  particular,  a  6  micrometer  thick  mylar  foil 

is  the  range  thickness  for  20  keV  electrons,  and  a  13  micrometer  Titanium  foil 

62 

is  the  range  thickness  for  50  keV  electrons.  The  experimental  setup  and 

6  2 

typical  results  for  both  foils  are  shown  in  Figure  3-4.  It  is  readily 
apparent  that  the  energetic  beam  pulse  full-width-half-maximum  is  approximately 
5  ns  and  that  the  output  current  peak  is  5.0  kA  with  a  6^  Mylar  foil  and 
4  ns,  0.64  kA  with  a  13p  Ti  foil.  These  presence  of  lower  currents  through  the 
Ti  foil  may  be  due  to  the  range  thickness  of  the  foil.  In  both  foils,  a  hole  is 
made  where  the  beam  passes  through,  indicating  much  of  the  beam  energy  is 
deposited  in  the  foils  raising  the  temperature  above  the  melting  point. 
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The  dramatic  difference  between  the  time  scales  of  the  Rogowski  coil 
current  diagnostic  and  the  Faraday  cup  diagnostic  indicates  the  presence  in  the 
Rogowski  coil  results  of  plasma  current  induced  in  the  plasma  channel  created  by 
the  beam.  In  order  to  investigate  this  phenomenon,  the  Faraday  cup  is  connected 
to  a  1.27  cm  diameter  aluminum  rod  and  inserted  into  a  drift  tube  of  10  cm  inner 
diameter  and  90  cm  length.  The  distance  of  the  Faraday  cup  from  the  center 
electrode  beam  port  is  adjustable  by  changing  the  length  of  the  rod.  The  plasma 
channel,  Al  rod,  and  drift  tube  may  be  considered  a  lumped  inductive  path.  The 
plasma  channel  has  resistivity.  Simple  Spitzer  resistivity  for  a  channel  of 
typical  parameters  is  given  by  the  equation: 

ri  -  0.0103  ln(A)  T(eV)‘^/^  -  0.013  fi-cm  .  (3-3) 

Assuming  typical  parameters,  ln(A)  -  10,^^  T  -  4  eV,^^  channel  cross  section  - 
0.95  cm  (estimated  from  a  witness  plate  exposed  1.25  cm  from  the  beam  port  in  a 
0.5  Torr  Ar  shot  shown  in  Figure  3-5)  for  a  channel  length  of  12.5  cm  the  net 
resistance  is  0.171  ohms.  Figure  3-6  is  an  open  shutter  photograph  of  the  beam 
induced  channel.  The  inductance  of  the  current  path  is  estimated  assuming  the 
current  remains  on  axis,  L  -  376  nH.  The  L/R  decay  time  of  this  lumped  circuit 
model  is  2.20  /iS .  This  value  is  comparable  to  the  experimental  values  found 
for  the  decay  constant  which  are  plotted  in  Figure  3-7.  This  is  an  indication 
of  the  relatively  large  cross  section  for  ionization  of  the  gas. 

ELECTRON  BEAM  CALORIMETER 

In  order  to  better  understand  the  electron  beam  production  in  a  plasma 

focus  device,  a  graphite  block  calorimeter  is  constructed . Other 
6  8 

researchers  have  reported  large  efficiencies  in  conversion  from  plasma  energy 
to  optical  energy.  Perhaps  even  though  the  beam  parameters  measured  at  the  end 
of  the  center  electrode  (beam  port)  indicate  relatively  low  e -beam- production 
efficiency  about  1  percent,  the  total  beam  energy  may  be  much  higher  but  lost  to 
the  walls  of  the  center  electrode.  Much  of  this  component  of  the  beam  may  be  of 
the  lower  energy  electrons  which  may  be  more  efficient  in  producing  optical 
radiation.  Hence,  we  wish  to  design  a  calorimeter  capable  of  measuring  the 
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total  beam  energy  at  the  entrance  to  the  center  electrode  as  well  as  along  its 
entire  length. 

A  calorimeter  is  constructed  (see  Figure  3-8)  from  a  graphite  block  to 
which  a  YSI44111  thermistor  is  thermally  connected.  The  thermistor  is  a 
negative  temperature  coefficient  resistor  which  is  capable  of  operating  in  the 
range  40°  to  150°  C.  Its  resistance  at  room  temperature  is  100  kohms .  The 
graphite  block  is  insulated  with  shrink  tubing  over  most  of  its  length  to 
provide  thermal  insulation  and  to  force  the  beam  current  to  pass  through  the 
entire  length  of  the  graphite.  Fingerstock  provides  electrical  contact  between 
the  block  and  ground. 

If  we  consider  only  the  graphite  (neglecting  heat  capacity  of  the  shrink 
tubing,  thermal  compound,  and  thermistor),  we  can  find  the  beam  energy  per 
change  in  temperature  calibration  factor.  Taking  the  graphite  mass  12.0  gm  and 
multiplying  by  the  heat  capacity  (Cp)  of  graphite,  0.712  J/gm-K,  we  obtain 
8.54  J/K.  Variations  of  Cp  with  temperature  are  negligible  in  the  range 
considered.  The  largest  error  in  the  experiment  may  be  the  loss  of  electron 
beam  around  the  outer  edge  of  the  calorimeter.  The  thermistor,  which  is  placed 
in  the  center  of  the  block,  provides  a  very  sensitive  method  of  measuring 
temperatures.  In  order  to  produce  a  record  of  these  changes,  the  thermistor  is 
connected  to  a  bridge  circuit.  This  circuit  is  shown  in  Figure  3-9.  A  Heath 
Schlumberger  Model  SR-255B  Chart  Recorder  is  connected  between  a  and  b.  The 
voltage  here  depends  on  the  change  in  the  resistance  of  the  thermistor: 


V  -  V  (  - = -  )  -  ( 

''ab  ^  ^  1+  R^^/R  ’  ^ 


Within  a  small  error: 


69 


(3-4) 


^1 

*4'2 


exp[/9^(— 


(3-5) 
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we  obtain: 


ab 


V  [ 


U  ^  -  H> 

K 


1+ 


Si 


(3-6) 


where  /9^  is  calculated  from  data  supplied  by  the  thermistor  manufacturer,  ■ 

4070  K.  Figure  3-10  shows  a  plot  of  output  voltage  versus  temperature 
difference.  Since  the  output  voltage  depends  on  the  intitial  balancing  of  the 
bridge  prior  to  a  shot,  there  is  a  weak  dependence  on  initial  temperature. 


In  order  to  interpret  the  raw  data  obtained  from  the  calorimeter 
diagnostic,  it  is  necessary  to  have  a  model  by  which  the  chart  recorder  results 
can  be  understood.  In  the  ideal  case,  a  short  pulse  of  power  P  -  IV  is  input  to 
the  calorimeter  as  in  Figure  3. 11(A).  The  duration  of  the  pulse  is  negligibly 
small.  The  calorimeter  diagnostic  is  time  integrating.  The  chart  recorder 
output,  which  is  proportional  to  total  energy,  should  be  a  step  function  as  in 
Figure  3.11(B).  In  actuality,  the  calorimeter- thermistor  system  does  not 
respond  instantaneously .  The  risetime  may  be  many  seconds.  The  exact  value 
depends  on  the  environment  immediately  surrounding  the  thermistor.  Calibration 
of  the  calorimeter  in  water  heat  baths  give  a  risetime  of  19s.  This  lies  within 
the  range  supplied  by  the  thermistor  manufacturer  (2.5s  in  well  stirred  oil, 

25s  in  still  air).  Taking  this  risetirae  into  consideration,  the  true  response 
to  a  pulse  input  may  be  expressed  as  AQ(l-e''/^r  )  where  Aq  is  proportional  to 
the  energy  and  tj,  is  the  risetime  of  the  thermistor.  See  Figures  3-ll(C).  The 
last  factor,  the  falltime,  is  due  to  heat  loss  from  the  calorimeter  through 
surrounding  gas,  shrink  tubing,  epoxy,  and  fingerstock.  The  heat  flow 
equation®^ is: 


_S_ 


S 

n 


-k 


dT(r,t) 
dn  ’ 


(3-7) 


where  q  is  the  heat  transfer  rate,  n  is  a  unit  vector  in  a  particular  direction, 
is  the  area  perpindicular  to  n,  T(r,t)  is  the  temperature,  and  k  is  the 
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thermal  conductivity  of  the  medium.  In  one  dimension,  we  can  integrate  from 
to  X2  to  obtain: 


q  -  kS 


(^1  -  ^^2^ 

(Xj^  -  x^)  ’ 


(3-8) 


where  2^'  cylindrical  coordinates,  the  corresponding  equation 

integrated  from  r^^  to  r2  over  a  z  length  Z  is; 


q  -  2»rkZ 


(T^  -  T2) 
ln(r2  /r^) 


(3-9) 


For  the  graphite  block,  the  heat  energy  stored  is  mg  Cp(T  -  where  mg 

Cp  is  the  product  of  mass  and  heat  capacity  of  the  graphite.  The  heat  energy 
flow  rate  from  the  calorimeter  is; 


q-  -  [m^Cp  (T-T_,,^^)]  . 


(3-10) 


Setting  this  equal  to  the  expression  above  and  assuming  room  temperature  as  the 
initial  temperature  bath,  we  arrive  at  the  equation: 


hT  -  ATq  e'^'^/^f^  , 


where  is  the  maximum( initial)  temperature  difference  and: 


(3-11) 


tj  -  m  Cp 


ln(r2  /r^) 
2»rkZ 


(3-12) 


The  quantity  t£  can  be  thought  of  as  a  falltime  that  must  be  incorporated  into 
the  calorimeter  response  which  may  be  described  by  the  equation; 


Aq(1  -  e'^^/*^r^  )  e'^^^^f^  . 


(3-13) 


See  Figure  3-ll(D).  Since  there  are  many  unknowns  in  the  calculation  of  t£, 
edge  corrections,  effective  fingerstock  conduction,  etc.,  we  will  only  use  the 
value  of  t£  obtained  experimentally  from  the  output  of  the  chart  recorder.  The 
quantity  Aq,  or  equivalently  ATq,  can  then  be  extrapolated  from  the  raw  data. 
The  total  energy  absorbed  by  the  graphite  is  m  Cp  AT^  where  AT^  near  room 
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temperature  is  0.183  K/mV  AV;  hence,  total  energy  per  output  voltage  calibration 
factor  is  1.56  J/mV. 

The  energy  absorbed  by  the  graphite  at  differing  total  inductive  energy 
during  focus  is  shown  in  Figure  3-12.  The  error  bars  are  representative  of 
uncertainty  in  t^  and  t£  measurements .  Notice  that  only  argon  gas ,  which  shows 
best  reproducibility  of  shots  as  judged  by  the  transmission  line  voltage  probe, 
shows  fairly  well  defined  slopes  at  the  two  filling  pressures  tested.  The  lines 
are  indicative  of  the  efficiency  of  inductive  energy  conversion  into  beam 
energy,  i.e.,  0.6  percent  at  2  Torr  and  1.4  percent  at  0.1  Torr  Ar.  The  most 
convenient  feature  of  this  diagnostic  is  that  it  may  be  inserted  into  the  hollow 
center  electrode  so  that  we  may  study  the  energy  transport  along  this  length. 
Figure  3-13  shows  how  the  beam  energy  is  transported  for  a  pressures  of  2  and 
0.5  Torr  argon.  The  solid  line  indicates  the  uniform  beam  divergence  line.  The 
lack  of  agreement  shows  the  importance  of  space  charge  effects  of  beam 
propagation  down  the  tube. 

ELECTRON  ENERGY  MEASUREMENTS 


A  simple  magnetic  analyzer,  based  on  principles  common  to  many  analyzers 

used  by  others,  has  been  designed  with  a  view  towards  compactness  and 
72-75 

versatility.  In  a  uniform  magnetic  field,  the  electron  will  follow  a 

circular  path  with  the  Larmor  radius  rj^  -  p/eB,  where  p  is  the  relativistic 
mechanical  momentum,  B  is  the  magnetic  field,  and  e  is  the  electronic  charge. 
The  kinetic  energy  T  corresponding  to  the  radius  r^^  is  given  by: 


T 


[(ecBr^^)^ 


+  (m 


-me 

e 


(3-14) 


where  c  is  the  velocity  of  light  and  m  is  the  electron  rest  mass.  Measurement 
of  the  Larmor  radius  is  accomplished  in  the  semicircular  focusing  electron 
energy  analyzer  shown  in  Figure  3-14. 


The  double  slit  collimation  system  in  Figure  3-14  limits  the  acceptance 
angle  and  width  of  the  electron  beam,  i.e.,  system  resolution.  The  thickness  of 
the  material  used  for  the  slits  in  the  collimator  should  be  greater  than  the 
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range  of  the  most  energetic  electrons.  For  example,  in  copper,  the  range  of 

2 

1.5  MeV  electrons  is  approximately  1  gm/cm  .  At  a  density  of  approximately 
9  gm/cm  ,  a  1.1  mm  thickness  of  copper  is  used. 

A  uniform  magnetic  field  of  0.53  Tesla  is  produced  by  a  single  rare  earth 
cobalt  magnet  of  dimension  2.54  x  2.54  x  0.64  cm  with  a  0.23  cm  separation  from 
a  symmetric  piece  of  soft  iron.  The  detectable  energy  range  is 
20  KeV-1.5  MeV.  The  detector  used  in  this  spectrometer  is  an  approximately  100 
micron  thick  coating  of  phosphor  P31  deposited  on  an  aluminum  substrate.  Since 
the  plasma  focus  is  operated  with  a  Torr  range  pressure  in  the  diode  region,  it 
is  desireable  to  provide  for  the  maintenance  of  vacuum  in  the  energy  analyzer  in 
order  to  avoid  an  electron  beam- filling  gas  interaction  which  could  interfere 
with  measurements.  A  fast  gate  valve  operated  by  a  solenoidal  plunger  system  is 
ideal  for  this  purpose.  See  Figure  3-15.  The  collimator  is  placed  immediately 
after  the  gate  valve.  Upon  the  opening  of  the  valve,  the  electron  beam  is 
injected  through  the  collimator  and  travels  a  semi-circular  trajectory  to  the 
detector.  The  results  of  each  shot  are  recorded  on  35  mm  black  and  white  film 
by  taking  an  open  shutter  photograph  of  the  phosphor  during  exposure  to  the 
electron  beam.  A  microdensitometer  scan  of  the  developed  film  yields 
information  about  the  energy  distribution  within  the  beam.  The  microdensito¬ 
meter  scans  the  optical  density  of  the  film,  which  may  not  necessarily  be  linear 
to  the  phosphor  exposure;  hence,  quantitative  information  about  the  energy 
spectrum  is  not  accurately  determined  by  this  method.  In  obtaining  the  energy 
data,  it  is  assumed  that  the  distance  from  the  center  of  the  collimator  slit  to 
a  point  on  the  phosphor  where  light  is  emitted  is  the  radius  of  the  electron 
orbit  in  the  uniform  field.  Such  an  assumption  is  rigorously  true  for  perfect 
collimation,  a  uniform  magnetic  field  configuration,  and  electron  orbits 
unaffected  by  space-charge.  Under  these  ideal  conditions,  we  may  convert  from 
distance  to  energy  measurement  by  writing  the  energy  in  convenient  laboratory 
units: 

T-  [(0.015Bs)^  +  (0.511)^]^'^2  .  q.SII  (MeV)  ,  (3-15) 

where  B  is  the  magnetic  field  in  kG  and  s  is  the  diameter  of  the  orbit  in  mm. 

The  resolution  of  this  spectrometer  is  discussed  in  Appendix  A. 
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Figures  3-16  and  3-17  show  results  for  a  typical  series  of  Ar  filling  gas 
shots.  Table  3-1  summarizes  peak  energies  obtained  for  argon  and  nitrogen  gas 
series  as  determined  by  the  sensitivity  of  the  detection  system,  i.e.,  amount  of 
phosphor  light,  f-stop  of  camera,  and  film  sensitivity.  These  were  done  on  a 
different  day  than  the  results  shown  in  Figures  3-16  and  3-17.  An  obvious 
notable  characteristic  is  the  irreproducibility  of  the  beam  energy  output  of  the 
device.  In  fact  it  was  observed  that  in  many  instances,  nitrogen,  gas  did  not 
show  any  result.  This  was  even  more  pronounced  in  the  cases  of  helium  and 
hydrogen. 


TABLE  3-1.  COMPARISON  OF  Ar  AND  N2  SERIES 


GAS  TORR 

PF  CURRENT (KA)  PEAK  ENERGY  OF  BEAM(KEV) 

Ar 

0.2 

55 

114 

Ar 

0.5 

127 

123 

Ar 

0.8 

144 

217 

Ar 

1.6 

137 

123 

Ar 

2.0 

137 

118 

N2 

0.3 

130 

85 

N2 

0.5 

137 

121 

CM 

1.0 

144 

none 

ELECTRON  BEAM  EMITTANCE 

Beam  quality  measurements  are  often  described  in  terms  of  emit-tance  which 
is  defined  as  the  area  occupied  by  the  beam  in  transverse  trace  space 
(coordinates  -  x ,y , dx/dz , dy/dz) ,  J  dx  dy  dx'  dy’  irregardless  of  the  density 
distribution  f^(x,y ,x' ,y' ) .  In  many  cases  the  actual  density  distribution  may 
be  defined  by  a  Gaussian  function  which  does  not  have  finite  limits;  hence,  the 
emittance  would  be  infinite.  In  addition,  emittance  is  an  invariant  quantity  in 
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systems  without  space-charge  interactions^^  (Liouville's  theorem);  therefore, 

erosion  of  beam  quality  due  to  f ilamentation  would  not  be  represented.  In  order 

to  acccv’.nt  for  these  asnects  of  beam  quality,  root -mean- square  (rms)  em' Ltar.ce 

75-78 

has  been  defined  as: 


t 

rms 


2  2 

4  [  <x  Xx '  > 


-  <xx 


(3-16) 


where  <4>>  represents  the  value  defined  by: 

<4>>  ”  ^  S  (x,y,x',y')  dx  dy  dx'  dy' 


(3-17) 


where  N  -  f^  dx  dy  dx'  dy' . 

79 

It  has  been  pointed  out  by  Lawson  that  the  rms  emittance  of  a  beam  is 
closely  related  to  its  entropy. 


Experimental  determination  of  f^(x,y,x' ,y' )  or  f2(x,x')  -  J  f^(x,y ,x' ,y' ) 
dy  dy'  may  be  accomplished  in  a  number  of  ways.  The  most  simplistic  may  be  two 
analyzers  each  having  an  aperture  through  which  a  portion  of  the  beam  may 
pass.  Emittance  meters  of  this  double  aperture  type  may  be  classified  into  one 
of  four  varieties  depending  on  the  types  of  apertures  in  the  two  analyzers, 
i.e.,  hole-hole,  hole-slit,  slit-hole,  or  slit-slit.  Let  us  consider,  for 
example,  a  hole-hole  emittance  meter.  The  upstream  analyzer  selects  a  certain 
x,y  position  within  the  beam  allowing  only  a  small  beamlet  to  pass  through  it. 
Then  as  the  downstream  analyzer  is  moved  about,  only  the  beamlet  particles  with 
slopes  x' ,y'  determined  by  the  relative  positions  of  the  two  analyzers  may  pass 
through.  The  intensity  distribution  f^  may  be  monitored  by  a  current  collector 
placed  immediately  behind  the  downstream  analyzer.  A  drawback  to  this  technique 
is  that  one  requires  a  continuous  beam  or  one  that  is  continuously  reproducible, 
over  many  shots,  which  is  often  not  the  case  in  the  plasma  focus  device.  In 
order  to  characterize  the  beam  in  a  plasma  focus  device,  we  require  an  emittance 
meter  capable  of  obtaining  single  shot  data.  A  single  shot  emittance  meter®*^’®^ 
has  been  designed  for  use  on  Febetron  and  an  IREB  device.  The  meter  used  here 
employs  a  radiachromic  film  detector  placed  4.32  mm  downstream  from  a  series 
of  400/im  wide  parallel  slits  placed  5.5  mm  apart.  See  Figure  3-18.  The  film 
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darkens  subsequent  to  exposure  to  the  electron  beam;  its  optical  density  is 
proportional  to  exposure  up  to  an  optical  density  of  approximately  1.0.  The 
film  may  be  rc'r!’'ned  with  an  optical  densitometer  to  determine  emittance.  This 
emittance  meter  may  be  classified  as  a  slit-hole  emittance  meter,  since  the  film 
is  scanned  through  the  pinhole  of  the  microdensitometer.  It  should  be  pointed 
out  that  of  the  above  mentioned  four  varieties  of  emittance  meters,  only  the 
slit-hole  meter  is  somewhat  less  straightforward  to  analyze.  In  order  to 
facilitate  analysis,  certain  assumptions  are  made  which  allow  us  to  require  only 
the  scan  along  the  central  y-axis  of  the  exposed  film.  The  measured  intensity 
pattern  is  proportional  to  J  f^(x,Ly' ,x' ,y' )  dy' .  If  f^  is  a  slowly  varying 
function  of  y  within  a  range  4y  -  +/-  Lo  where  a  is  the  rms  value  of  y' ,  then 
the  measured  intensity  pattern  is  approximately  equal  to  the  one  obtained  by  the 
hole-slit  meter, 

f2(x,y,x')  -  J  f^(x,y,x' ,y' )  dy'  .  (3-18) 

This  condition  is  easily  met  in  beams  with  a  small  transverse  temperature  or 

meters  with  a  small  slit-detector  distance,  L.  In  addition,  our  analysis  will 

consider  only  axisymmetric  beams,  i.e.,  all  measurable  quantities  are  functions 
2  2 

of  r  -  y(x  +  y  )  only.  Also,  zero  beam  rotation  is  assumed  with  a  thermal 
(Maxwellian)  distribution  in  x*  -  dx/dz  and  y'  -  dy/dz  space.  These  may  be 
quite  reasonable  assumptions  for  a  real  beam  in  the  absence  of  external  magnetic 
field.  Under  these  assumptions,  the  overall  trace  space  distribution  function 
f^  may  be  given  by: 

f^(x,y,x' ,y' )  -  g(x,y)  exp-(( (x' -x' )^+(y' -y' )^]/2a^(x,y)]’  ,  (3-19) 

here  g(x,y)  -  g(r)  is  the  x-y  spatial  dependence  of  the  distribution  function, 
‘^(x.y)  ”  o(r)  is  the  rms  value  of  the  angle  (velocity)  in  x'-y'  space  given  by: 

y{/[(x'-x')^+(y'-y')^]f^dx'dyV/f^dx'dy']’  ,  (3-20) 

and  x'(r)  and  y' (r)  are  the  mean  angles  in  x'-y'  space  given  by: 

J  (x'  or  y')  f^  dx'  dy'  /  f  dx'  dy'  ,  (3-21) 
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which  are  typically  due  to  overall  beam  divergence  or  convergence.  We  have 
written  x,y  expli''itly  in  the  above  equations  to  point  out  the  fact  that  the 
film  is  scanned  along  the  y-0  axis  (actually,  we  assiome  the  film's  y-0  axis  is 
the  same  as  the  upstream  analyzer's  y-0  axis)  and  cylindrical  symmetry 
determines  the  parameters  g(x,y),  a(x,y),  x'(x,y),  and  y'(x,y)  in  full  x-y 
space . 


An  analytical  method  of  analyzing  experimental  data  produced  by  a  slit-hole 

fi3 

emittance  meter  shall  be  presented.  The  slits  at  Xj^  on  the  upstream  analyzer 
produce  sheet  beamlets  which  impact  on  the  detector  revealing  their  individual 
angular  distributions  by  the  amount  of  spread  as  shown  in  Figure  3-19.  This  is 
measured  quantitatively  with  the  microdensitometer  which  shows  typically  a 
gaussian  profile  with  peak  height  rms  width  and  displacement  xj^. 
Continuous  functions,  ^(x),  a(x) ,  and  x'(x),  may  be  constructed  from  the 
discrete  sets  of  data.  In  fact,  a  computer  program  has  been  written  which  will 
fit  fi/a  to  a  gaussian,  a  to  a  second  order  polynomial,  and  a  -  x'(x)  r/x  to  a 
linear  function.  Since  all  quantities  are  assumed  to  be  measured  along  y-0,  we 
find  due  to  symmetry  that  /9(r)-^(x),  a(r)-a(x),  and  Q(r)-a(x).  It  has  been 

O  0 

shown  that  by  utilizing  a  transformation  of  the  type: 


ycr  -X  ) 

the  quantities  needed  for  emittance  may  be  found  as: 


<x2>  _  ^  2|-c_^ar_^ 

N  Jo  Jx  2  2,  ’ 

y(r  -X  ) 


N  Jo^Jx  ,.2  2,  Jx  2  2  J  ■ 

y(r  -x  )  rycr  -x  ) 


<xx'>  -  ^  r  r  ^  dx 

^  N  Jo  Jx  2  2,  ' 

y(r  -X  ) 


where , 


(3-22) 


(3-23) 


(3-24) 


(3-25) 


3-12 


NSWC  TR  88-394 


^  ^  -C- 


Bar  dr 

//  2  2, 

y(r  -X  ) 


dx 


(3-26) 


Figure  3-20  shows  a  microdensl tometer  sc^n  of  a  radiachromic  film  which  was 
exposed  at  the  end  of  the  center  electrode  as  it  leads  into  the  10  cm  diameter 
drift  tube.  The  shot  was  a  750  mTorr  Ar  gas  shot  with  a  plasma  current  of  150 
kA  at  focus.  For  this  shot,  the  radiachromic  film  was  covered  with  a 
0.00025 -inch  thick  aluminum  foil  to  prevent  ultraviolet  or  low  energy  electrons 
from  exposing  the  film.  Only  three  slits  across  the  meter  were  exposed  since 
the  hollow  center  electrode  diameter  is  14.4  mm  and  the  slit  separation  is  5.5 
mm.  It  should  be  emphasized  here  that  this  measurement  was  performed  at  the 
center  electrode  beam  port  since  the  beam  emittance  may  vary  as  the  beam 
propagates  further  downstream  due  to  scattering,  space  charge,  etc.  The  a  value 
at  the  center  of  the  beam  is  630  ~rad  and  at  x  -  5.5  mm,  the  value  is  550 
mrad.  The  value  of  a  fit  to  a  straight  line  is  0.021r,  and  the  overall  spatial 
distribution  is  approximately  a  flat  top.  From  Equations  (3-23)  through  (3-25) 
we  find  the  rms  emittance  for  a  flat  top  spatial  distribution  with  a  Gaussian 
transverse  velocity  profile  to  be  given  by: 


2boQ  ,  (3-27) 

where  b  is  the  radius  of  the  flat  top  distribution  and  a^  is  the  peak  of  the  rms 
width  function  a.  Plugging  in  the  values  b-7.2  mm,  560  mrad,  we  get  8060 
mm-mrad.  A  computer  program  has  been  written  (Won  Namkung)  which  enables  us  to 
numerically  Integrate  the  equations  for  the  rms  values.  The  program  results 
show  a  rms  emittance  of  8040  mm-mrad.  A  two-dimensional  contour  plot  of 
f2(x,x')  is  shown  in  Figure  3-21. 
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FIGURE  3-1.  E-BEAM  ROGOWSKI  COIL  TRACES  (INTEGRATED)  1.42  kA/div 
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FIGURE  3-1.  (CONT.)  E-BEAM  ROGOWSKI  COIL  TRACES  (INTEGRATED) 
i.42  kA/div. 
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FIGURE  3-3.  PLOT  OF  BANK  CURRENT  AT  FOCUS  I  ^  VERSUS  PRESSURE  FOR  A 
VARIETY  OF  FILLING  GAS 
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FIGURE  3-4. 
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FIGURE  3- 


FIGURE  3- 


O  ‘ 

c. 


WITNESS  FILM  FOR  ELECTRON  BEAM  (FIVE  EXPOSURES  ON  A  WITNESS 
FILM  FOR  E-BEAM  IN  0.5  TORR  Ar) 


OPEN  SHUTTER  PHOTO  OF  E-BEAM.  OPEN  SHUTTER  PHOTOGRAPH  OF 
BEAM  PASSING  THROUGH  0.32  TORR  Ar.  IN  THE  CENTER  OF  THE 
PHOTO  IS  A  GRAPHITE  BLOCK  CONNECTED  TO  1/2 -INCH  Al  ROD 
FOR  FARADAY  CUP  MEASUREMENTS. 
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FIGURE  3-7.  PLOT  OF  ELECTRON  BEAM  INDUCED  PLASMA  CURRENT  DECAY  TIME  AS  A 
FUNCTION  OF  GAS  PRESSURE  AND  DISTANCE  WHERE  THE  FARADAY  CUP 
IS  12.5  AND  33  cm  FROM  THE  BEAM  PORT 
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FIGURE  3-8.  E-BEAM  CALORIMETER 


PLOT  OF  OUTPUT  VOLTAGE  VS  TEMPERATURE  CHANGE 
FOR  THE  CALORIMETER 
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FIGURE  3-10.  OUTPUT  VOLTAGE  OF  BRIDGE  AS  A  FUNCTION  OF  THERMISTOR  CHANGE 
IN  TEMPERATURE  FOR  THE  APPLIED  VOLTAGE  5.05  V 
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FIGURE  3-11.  ANALYSIS  OF  CALORIMETER  RESPONSE 
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HOLLOW  CENTER  ELECTRODE 
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FIGURE  3-14.  MAGNETIC  ELECTRON  ENERGY  ANALYZER 
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FIGURE  3-16.  MICRODENSITOMETER  SCAN  OF  PHOTOGRAPHIC  FILM  SHOWING 
THE  ELECTRON  ENERGY  DISTRIBUTION 
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FIGURE  3-17.  MICRODENSITOMETER  SCAN  OF  PHOTOGRAPHIC  FILM  FOR  ENERGY 
DISTRIBUTION  OF  A  SERIES  OF  Ar  SHOTS 
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RADIACHROMIC  FILM 
/ 


16  X  195  mil  Ta  BARS 
20  mil  SLIT  WIDTH 


FIGURE  3-18.  CROSS-SECTIONAL  VIEW  OF  THE  SLIT- PINHOLE  EMITTANCE  METER 
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FIGURE  3-19.  A  TYPICAL  GAUSSIAN  INTENSITY  PROFILE  CREATED  BY  THE 

ANALYZER  SLIT  AT  Xj.  FOR  L-1 ,  THE  MEAN  ilVERGENCE  ANGLE 

X'  rms  WIDTH  a.,  AND  HEIGHT  /9.  ARE  SHOWN, 
i  1  1 
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MICRODENSITOMETER  SCAN  OF 
EMITTANCE  METER  RESULT 
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FIGURE  3-20.  MICRODENSITOMETER  SCAN  OF  RADIACHROMIC  FILM  FOR  EMITTANCE 
ANALYSIS 
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CONTOURS  OF  CONSTANT  f2(x,  x' ) 
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FIGURE  3-21.  fjCx.x')  CONTOUR  PLOT  SHOWING  0.01,  0.1,  0.5,  AND 
0.9  CONTOURS 
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CHAPTER  4 
ION  BEAM  STUDIES 


ION  BEAM  FLUX  MEASUREMENTS 

The  ion  beam  flux  is  measured®^  using  CR-39  nuclear  track  detector.®^ 

CR-39  plastic  is  damaged  by  the  passage  of  energetic  ions  as  they  enter  the 
plastic.  The  damaged  area  is  then  typically  etched  out  of  the  plastic  leaving  a 

QC 

visible  track.  The  etching  of  the  track  progresses  as  shown  in  Figure  4-1. 

The  damaged  area  etches  at  a  velocity  v^,  which  is  more  rapid  than  the  etch  rate 
of  the  undamaged  plastic,  v  The  sensitivity  of  this  detector  is  quite  high. 

O 

Protons  of  energies  down  to  several  keV  have  produced  tracks  in  etched  CR-39, 
Lighter  ions,  such  as  protons,  tend  to  produce  tracks  with  a  very  small  damage 
radius  whereas  heavier  ions  will  produce  a  larger  radius  of  damage  due  to  their 
higher  atomic  number,  but  a  shallower  track  due  to  their  shorter  range  in  the 
plastic.  This  phenomenon  is  observed  quite  frequently  in  Thomson  spectrograms 
which  will  be  discussed  later  in  this  chapter. 

In  order  to  measure  the  flux,  a  copper  mask  is  made  with  small  holes  along 
a  radius  every  7  mm.  The  CR-39  detector  is  placed  immediately  behind  this  mask 
and  exposed  to  the  plasma  focus  discharge.  The  CR-39  may  be  rotated  enabling  up 
to  12  exposures  without  opening  the  system  to  vacuum.  After  the  CR-39  is 
exposed  to  the  focus,  it  is  etched  in  6.25  N  solution  of  NaOH  at  70°C  for  30 
minutes.  At  this  time,  the  damage  caused  by  the  passage  of  the  ions  through  the 
plastic  has  been  etched  sufficiently  to  be  barely  visible.  The  resulting  track 
density  is  much  too  high  to  enable  a  measurement  of  density  under  a  light 
microscope.  A  Formvar  replica  of  the  track  plate  is  taken  so  that  scanning 
electron  microscopy  may  be  used.  The  Formvar  is  in  a  2  percent  solution  of 
ethelyne  dichloride  so  that  it  may  be  dropped  onto  the  track  plate  with  a 
dropper.  The  evaporation  of  the  solvent  leaves  only  a  thin  layer  of  Formvar 
covering  the  CR-39  and  filling  the  etched  track  pits.  When  this  layer  is  peeled 
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from  the  CR-39,  a  negative  image  of  the  tracks  has  been  formed.  These  track 
spikes  are  much  easier  to  observe  under  the  scanning  electron  microscope  than 
the  corresponding  track  pits  and  enable  us  to  count  the  track  density  very 
simply. 

The  number  of  ions  per  unit  solid  angle  is  obtained  by  counting  the  tracks 
per  unit  area  on  the  picture,  nQ,  and  putting  into  the  equation: 

^  -  n^  (L^+R^)  /cosd  cosa^  ,  (4-1) 

where  M  -  picture  magnification,  L  -  distance  from  focus  to  pinhole  plate,  R  - 

radius  from  axis  of  system  to  pinhole,  6  —  tan  ^(r/L),  and  a  ,  —  tilt  angle  of 

dN  ^ 

sample  on  scanning  electron  microscope  stage.  The  flux,  —  is  found  from  the 

above  equation  and  plotted  as  a  function  of  polar  angle  in  Figure  4-2  for  a  250 

mTorr  N2  shot  and  a  250  mTorr  H2  shot.  On  this  shot  we  differentiated  between 

proton  tracks  and  heavier  ion  impurity  tracks  such  as  nitrogen  by  the  shape  of 

the  etched  track  pit.  Impurities  in  the  focus  are  important  constituents  of  the 

ion  beam  especially  at  larger  angles.  This  is  consistent  with  the  hypothesis 

that  the  heavier  impurities  because  of  their  greater  mass/charge  ratio  will  be 

on  the  outside  of  the  sheath  as  it  collapses.  Because  of  larger  larmor  orbits, 

/  0 

these  ions  may  be  preferentially  accelerated  at  larger  polar  angles. 

Integration  of  the  dN/dO  may  be  approximated  by  a  summation  (assuming  axial 
symmetry)  to  estimate  the  total  number  of  ions  accelerated  within  the  range 
tested. 


X  2x  ^  sine  Ae  ,  (4-2) 

gives  1.05  X  10^^  particles  within  a  22-degree  cone  for  the  nitrogen  exposure 
shown  in  Figure  4-2. 

We  attempted  to  see  some  filling  gas  dependence  by  finding  the  track 
density  on  axis  at  a  distance  of  8.3  cm  from  the  tip  of  the  center  electrode. 

250  mTorr  hydrogen,  500  mTorr  helium,  250  mTorr  nitrogen,  and  500  mTorr  Ar  shots 
produced  1.26,  1.05,  5.39,  and  1.04  X  10^  tracks  per  cm^ .  The  plasma  currents 


4-2 


NSWC  TR  88-394 


for  these  shots  at  focus  were  135,  140,  145,  and  150  kA,  respectively.  It  is 
apparent  that  any  pressure  dependence  may  be  overwhelmed  by  shot  to  shot 
variation  of  the  focus.  In  addition,  hydrogen  results  at  pressures  60,  125, 

250,  500,  1000,  and  2000  mTorr  were  1.45,  2.56,  1.26,  1.80,  0.54,  and  2.14 
X  10^  tracks/cm^ . 

ION  PINHOLE  CAMERA  ANALYSIS 

89 

A  stereoscopic  pinhole  camera  is  constructed  in  order  to  observe  the 
spatial  distribution  of  the  ion  source  of  the  plasma  focus  device.  A  very 
important  consideration  in  the  design  of  such  a  camera  is  the  resolution  of  the 
system.  The  quantities  important  for  resolution  analysis  are  D  -  pinhole 
diameter,  0  -  object  to  pinhole  distance,  and  L  -  image  to  pinhole  distance. 

With  these  definitions,  it  is  a  simple  geometrical  problem  to  find  R^  the 
pinhole  camera  resolution  spot  size,  i.e.,  the  distance  between  two  objects 
which  can  just  be  resolved; 

-  D  (0/L  +1)  ,  (4-3) 

See  Figure  4-3. 

The  pinhole  camera  is  constructed  from  two  pinholes  of  0.075  mm  diameter  on 
a  25^m  thick  copper  foil  separated  by  10  mm  (-P) .  A  CR-39  detector  is  placed  15 
mm  downstream  of  the  pinholes.  When  the  camera  is  placed  75  mm  from  the  focus, 
Rj,  is  0.3  mm.  Figure  4-4  shows  several  pinhole  pictures  taken  of  the  plasma 
focus  at  250,  250,  and  500  mTorr  of  Ar  filling  gas  (top  to  bottom).  Left  and 
right  exposures  correspond  to  the  same  shot,  and  the  figure  shows  the  CR-39  as 
viewed  from  the  exposed  side. 

Assuming  straight  trajectories  from  the  source,  stereoscopic  'analysis  may 
be  applied  to  such  a  double  pinhole  camera.  The  geometry  is  shown  in 
Figure  4-5.  For  simplicity,  suppose  two  point  ion  sources  shown  at  different  x 
and  z  coordinate  emit  uniformly  over  the  angle  subtended  by  the  pinholes.  The 
x-axis  is  defined  by  the  line  passing  through  both  pinholes  and  the  z-axis  is 
defined  by  the  line  from  the  source  to  the  pinhole.  The  separation  between  the 
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two  images  will,  in  general,  be  different  at  the  left  and  right  images,  their 
separations  being  and  Sj^,  respectively.  It  is  straightforward  to  find  the 
separation  of  the  sources  in  z  and  x  coordinates,  Az  and  Ax, 

Az  -  O^(Sj^-S^)  /LP  ,  (4-4) 

Ax  -  0(Sj^+Sj^)  /2L  .  (4-5) 

We  may  identify  some  structures  common  to  both  left  and  right  pictures  in 
Figure  4-4  as  well  as  their  relative  locations.  The  difficulties  in  this 
analysis  are  the  lack  of  sharp  detail  in  most  of  the  structures  identified, 
which  typically  leads  to  an  error  of  +/-  several  millimeters,  hence  this 
stereoscopic  technique  may  be  of  limited  use  in  determining  the  three- 
dimensional  ion  source;  however,  the  two-dimensional  information  is  quite 
reliable . 

Some  interesting  features  of  the  pictures  are  the  up  to  four  concentric 
circle  structures  at  effective  source  radii  of  1.47,  0.92,  0.60,  and  0.24  cm. 

In  addition,  at  least  one  point  source  appears  within  the  smallest  radius 
circle.  Similar  structure  has  been  observed  at  other  laboratories^®  with  no 
apparent  explanation.  It  must  be  remembered  that  this  is  a  time  integrated 
diagnostic.  The  point- like  source  at  the  center  may  be  comparable  in  size  to  the 
final  diameter  (sub-millimeter)  of  the  pinch  observed  in  streak  or  framing 
photographs.  The  larger  concentric  circles  may  be  indicative  of  earlier  ion 
acceleration  as  the  cylindrically  symmetric  sheath  is  collapsing  toward  the 
axis . 

THOMSON  SPECTROMETER  ION  DIAGNOSTIC 

Q1 

A  compact  Thomson  spectrometer  designed  by  M.  J.  Rhee  is  used  to  analyze 
the  ion  beam  for  charge  state  resolved  energy  spectra.  The  Thomson  spectrometer 
is  a  charge- to -mass  ratio  and  energy  analyzer  based  on  particle  deflection  in 
parallel  electric  and  magnetic  fields.  Reference  91  is  a  complete  description 
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of  the  spectrometer  and  particle  deflection  analysis.  A  short  synopsis  will  be 
given  in  this  work  for  the  sake  of  completeness. 


The  spectrometer  is  based  on  two  high  induction  permanent  magnets.  A 
copper  foil  covering  the  magnet  surface  allows  it  to  be  biased  to  a  high 
voltage.  See  Figure  4-6.  In  the  electric  field,  charged  particles  will  follow 
a  parabolic  trajectory  with  a  net  deflection  angle: 


(4-6) 


where  Ze  is  the  net  charge,  jEdl  is  the  path  integral  of  the  electric  field 
along  a  particle  trajectory,  and  T  is  the  energy.  In  the  magnetic  field,  they 
will  follow  a  circular  trajectory  with  a  net  deflection  angle  which  is 
perpendicular  to  the  electric  deflection: 


(4-7) 


where  p  is  the  ion  momentum,  Jfidl  is  the  integral  of  the  magnetic  field,  and  Au 
is  the  product  of  ion  mass  number  and  unit  nucleon  mass.  The  resultant  pattern 
recorded  on  CR-39  nuclear  track  detector  is  a  set  of  parabolas  corresponding  to 
each  charge- to -mass  ratio  swept  out  by  energy  spread  in  the  ion  beam  described 
by  the  well  known  parabola  equation: 


m 


Ze  (jBdl)^^ 


Au  lEdl 


(4-8) 


The  collimation  for  the  system  is  provided  by  two  pinholes.  The  first 
upstream  pinhole  with  1.5  mm  diameter  (-Dj^)  is  placed  50  cm  (-Lj^)  upstream  of 
the  spectrometer  in  a  fast  (ms)  gate  valve  which  maintains  vacuum  in  the 
spectrometer  region.  The  second  pinhole  of  lO/im  micrometer  diameter  (-D2)  is 
just  downstream  of  the  field  region.  The  CR-39  detector  is  placed  5  mm  (“L2) 
downstream  of  tne  second  pinhole.  The  resolution  of  such  a  system,  and  Thomson 
spectrometers  in  general,  is  described  in  Appendix  B. 
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The  calibration  of  the  magnetic  field  integral  is  done  with  an  alpha 

particle  source.  The  collimated  beam  is  passed  through  a  collimation  systemwith 
and  without  the  magnetic  field.  The  value  of  /fidl  is  found  to  be  6.06  X  10"^  T- 
m.  Direct  measurements  with  a  Hall  probe  agree  to  within  1.5  percent.  The 
integral  of  electric  field  is  found  by  reduction  from  Thomson  spectrograms  of 
known  ion  species,  jEdl  -  1.25  EL,  where  EL  is  the  electric  field  times  the 
length  for  the  ideal  parallel  plate  capacitor. 

A  typical  spectrogram  is  shown  in  Figure  4-7.  For  this  spectrogram  jEdl  - 

22.4  kV.  Some  interesting  features  to  note  are  the  reference  B-  axis  which  is 

made  by  exposing  the  spectrometer  without  applying  an  electric  field.  In 
addition,  there  always  exists  a  spot  at  the  origin  of  the  Thomson  spectrometer 
coordinates.  This  is  due  to  the  undeflected  passage  of  particles  through  both 
pinholes  as  would  be  expected  of  neutrals. 

Also  note  that  due  to  the  polepieces,  there  is  a  finite  acceptance  angle 

which  is  determined  by  the  aspect  ratio  of  the  gap,  0.270  rad.  The 

identification  of  ion  species  depends  on  the  assumption  that  all  charge  states 
below  the  one  in  question  exist  and  extend  to  a  constant  velocity  line.  The 
charge  state  is  found  simply  by  counting  the  equally  space  intersections  of 
parabolas  with  this  constant  velocity  line.  The  mass  is  found  from  known  charge 
and  charge-to-mass  ratio. 

The  energy  spectrum  may  be  obtained  by  directly  counting  the  particle 
tracks  in  CR-39.  Each  parabola  may  be  divided  into  small  segments  of 
constant  electric  deflection  angle  which  correspond  to  a  certain  AE^.  The 
number  of  tracks  in  each  segment,  AN^^/A^^  ,  are  counted  and  the  spectrum  is 
found  from  the  equation: 

AN.  A«. 

dN  1  1 

—  .  — -  — -  (U-9) 

dE  A^.  AE.  ^  ^ 

1  1 

The  energy  spectra  of  the  ion  species  obtained  are  shown  in  Figures  4-8 
through  4-12.  The  energy  uncertainty  is  discussed  in  Appendix  B. 
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One  major  disadvantage  of  the  Thomson  spectrometer  system  is  that  the 
collimation  system  by  its  double  pinhole  nature  selects  only  a  small  solid 
angleof  the  ion  beam.  For  example,  one  could  imagine  the  second  downstream 
pinholeselects  only  a  portion  of  the  image  formed  by  the  upstream  pinhole. 

Howe  ’r  in  all  cases  measured,  we  have  taken  the  best  spectrogram  from  a  sample 
of  10  or  more.  Another  disadvantage  of  the  Thomson  spectrometer  is  that  it  is 
not  sensitive  to  the  high-energy  tail  observed  at  many  laboratories. 

NEUTRAL  PARTICLES 

92 

The  presence  of  neutral  particles  is  investigated.  This  study  was 

O 

motivated  by  the  observation  that  more  than  10  particles  travel  through  the 

hollow  center  electrode  in  the  electron  beam  geometry.  A  piece  of  CR-39  was 

placed  several  cm  from  the  beam  port,  and  while  the  electron  beam  proved  very 

unstable  and  often  missed  the  piece  of  CR-39,  there  was  a  consistent  production 

of  tracks  on  the  CR-39  in  a  sharp  disc  as  if  the  beam  was  well  collimated.  As  a 

matter  of  fact,  the  tracks  behaved  as  if  the  beam  was  collimated  by  the 

acceptance  of  the  hollow  center  electrode.  Preliminary  evidence  by  magnetic 

deflection  indicates  that  the  ions  have  a  greater  proportion  of  neutrals  than 

9  3 

the  forvard  ions,  a  phenomena  obseirved  by  the  Stevens  group  also.  Typically 
in  the  forward  direction  there  are  20-30  percent  neutrals. 

EMITTANCE  OF  ION  BEAM 

A  slit-pinhole  emittance  meter  is  used  to  determine  the  rms  emittance  of 
the  plasma  focus  produced  ion  beam.  The  meter  is  similar  to  the  one  used  for 
the  electron  beam  emittance.  In  this  case,  six  slits  are  20^jm  in  width, 
constructed  from  stainless  razor  edge  blades.  The  separation  between  slits  is 
9.50  mm  +/-  0.05  mm.  A  CR-39  nuclear  track  detector  is  placed  8.36  mm 
downstream  of  the  slits.  The  emittance  meter  is  placed  50  cm  downstream  of  the 
focus.  Due  to  the  wide  ion  beam  dispersion,  we  define  a  useful  beam  solid  angle 
by  placing  a  1.74  cm  diameter  diaphragm  22.5  cm  downstream  of  the  focus  which 
defines  a  solid  angle  /2xsin^dfi  -  0.0047  srad  (or  a  2.2-degree  cone).  Sub¬ 
sequent  to  exposure  of  the  meter  to  the  beam  produced  by  a  1  Torr  Ar  gas  shot 
with  plasma  current  of  190  kA,  the  CR-39  is  etched  for  1  hour  in  6.25  N  NaOH  at 
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70°C.  The  sample  is  then  placed  on  a  microscope  stage  with  an  attatched  dial 
gauge  for  measurements  to  0.01  mm  precision.  A  50x50  sq.  /im  area  on  the 
trackplate  defines  the  pinhole  size.  The  number  of  particle  tracks  within  this 
area(pinhole)  are  counted  and  plotted  as  a  function  of  tsi  in  Figure  4-13. 

These  six  individual  peaks  are  fit  to  a  Gaussian  function  to  find  their  mean 
widths.  It  was  found  that  to  within  a  very  slight  variation,  the  widths  {o)  are 
constant  as  would  be  expected  from  a  cold  beam  with  a  source  of  finite  spatial 
extent.  The  peak  values  of  the  distribution  are  fit  to  a  Gaussian  function  to 
find  the  mean  overall  beam  width,  and  the  beam  divergence  is  found  as  explained 
in  Chapter  3.  The  empirical  values  are  -  69.7  mrad,  b  -  19.9  mm,  and  a  - 
r/0.217.  Figure  4-14  is  a  two-dimensional  contour  plot  of  f 2 .  The  rms 
emittance  is  calculated  from  Equations  (3-23)  through  (3-25)  and  found  to  be 
4baQ  for  a  distribution  described  by  the  equations: 

2  2 

a(r)  -  r/a,  ^(r)  -  ,  o(r)  -  .  (4-10) 

Hence  the  rms  emittance  is  found  to  be  5546  mm-mrad. 
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ORIGINAL  SURFACE 


/ 


FIGURE  4-1.  TRACK  PIT  ETCHING.  PATHS  OF  ION  THROUGH  THE  PLASTIC  WILL 
ETCH  AT  A  LARGER  VELOCITY.  V^.  THAN  BULK  ETCH,  v^.' 
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ION  FLUX  AS  A  FUNCTION  OF  POLAR  ANGLE 
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FIGURE  4-2.  ION  FLUX  DISTRIBUTION  AS  A  FUNCTION  OF  POLAR  ANGLE 
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figure  4-3. 


resolution  of  thf  t'TMii/-. 
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FIGURE  4-4.  PINHOLE  PICTURES  OF  250.  250,  AND  500  mTORR  Ar  SHOTS 
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TOP  VIEW 


0  5  cm 


FIGURE  4-6.  THE  COMPACT  THOMSON  SPECTROMETER 
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FIGURE  4- 


THOMSON  SPECTROGRAM  OF 
ARGON  SHOT 


ELECTRIC  DEFLECTION 
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Omrad  100  mrad 


AN  EXAMPLE  OF  A  TYPICAL  Ar  RESULT  FROM  THE  SPECTROMETER 
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4-16 
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FIGURE  4-9.  ENERGY  SPECTRA  FOR  He  IONS 
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FIGURE  4-10,  ENERGY  SPECTRA  FOR  N  IONS 
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FIGURE  4-11.  ENERGY  SPECTRA  FOR  Ne  IONS 
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ENERGY/CHARGE(keV/Z) 


FIGURE  4-12.  ENERGY  SPECTRA  FOR  Ar  IONS 


EMITTANCE  DATA 


FIGURE  4-13.  ION  TRACK  DISTRIBUTIONS  FOR  EMITTANCE  ANALYSIS 
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CHAPTER  5 

ANALYSIS  OF  EXPERIMENTAL  RESULTS 


NONREPRODUCIBILITY 

The  nonreproducibility  of  the  plasma  focus  device  is  a  common  problem  with 
most  pulsed  power  systems.  The  performance  of  the  plasma  diode  depends  on  the 
performance  of  a  plasma  instability  which  may,  by  its  nature,  be  very  sensitive 
to  system  parameters  such  as:  impurity  gases  which  have  leaked  into  the  system; 
electrode  contamination  with  pump  oil;  asymmetries  in  the  electrodes  and/or 
insulators;  etc.,  which  may  vary  randomly  on  a  shot-to-shot  basis.  These 
considerations  constrained  us  to  use  single  shot  diagnostics  which  showed  a  wide 
variation  of  results. 

The  electron  beam  current  measured  by  the  Rogowski  coil  shewed  wide 
variations  of  approximately  2  kA  at  lower  pressures  (less  than  0.2  Torr)  of  the 
filling  gas.  At  the  higher  pressures  the  results  were  more  consistently 
lower.  The  Faraday  cup  measurements  after  the  foil  were  somewhat  erratic  with 
many  shots  (60  percent)  showing  null  result.  This  was  apparently  due  to  the 
poor  beam  propagation  in  the  short  2  cm  length  of  10  cm  diameter  drift  tube 
before  the  foil.  In  fact,  in  many  experiments  the  beam  damaged  the  brass  plate 
which  separated  the  electron  beam  gun  from  the  10  cm  tube  indicating  the  beam 
made  a  180-degree  turn.  It  may  be  of  use  to  note  that  the  steady  state  space 
charge  limiting  current  for  a  100  keV  electron  beam  in  an  infinite  tube  of  this 
diameter  is: 

17  (7^/^-!)^'^^  /  l(l+21n(ryr^))(l-f  )]  (kA)  .  (5-1) 

Where  rjj  is  the  beam  radius,  r^  is  the  drift  tube  radius  and  f^  is  the 
fractional  charge  neutralization.  This  value  is  less  than  0.2  kA  for  our 
geometry  and  beam  energy  when  fg  is  0.  The  electron  spectrometer  results  showed 
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at  least  as  wide  a  variation  since,  in  addition  to  the  short  length  of  10  cm 
tube,  only  a  small  fraction  of  the  electron  beam  was  selected  by  a  1/16 -inch 
diameter  pinhole.  Many  shots  showed  no  phosphor  light  while  some  showed  quite 
good  results  for  the  same  parameters.  In  addition,  the  collimator  slits  became 
clogged  with  gr  iphite  over  a  period  of  time  which  required  cleaning. 

The  Thomson  spectrometer  diagnostic  also  showed  a  wide  range  of  results. 
This  was  explained  in  Chapter  4  as  being  due  to  the  nature  of  the  two  slit 
collimation  system  which  selects  out  only  a  small  portion  of  the  ion  beam.  It 
should  also  be  mentioned  that  there  were  series  of  shots  that  showed  no  current 
drop  or  voltage  signal  with  no  explanation  readily  apparent.  This  would  quite 
often  fix  itself  the  next  day  or  when  a  new  insulator  was  placed  in  the 
device.  This  was  characteristic  of  the  ion  geometry  more  often  than  the 
electron  gun  geometry.  This  behavior  may  be  attributed  to  the  conditioning  of 
the  insulator  or  electrodes  as  discussed  in  Chapter  2  and  is  different  from  the 
shot  to  shot  variation. 

ANALYSIS 

From  a  nitrogen  parabola  and  the  ion  flux  measurement,  an  estimate  of  the 
total  ion  current  can  be  made.  In  the  parabola  analyzed,  there  are  1436  N"*", 

870  N  ,  and  216  N  ions.  This  gives  an  average  charge  of  1.52  neglecting 
neutrals.  It  should  be  emphasized  that  this  may  lead  to  an  underestimate  of  the 
total  current  since  ions  have  undergone  charge  exchange  or  recombination  as  the 
pass  through  a  dense  plasma  and  a  neutral  gas  after  acceleration.  Current  is 
approximated  by  (total  number  of  ions)  times  (average  charge  per  ion)  times  e 
divided  by  (the  pulse  width).  If  the  pulse  width  is  estimated  to  be  the  same  as 
the  electron  beam  width  4  ns,  the  net  ion  current  is  6.4  amps,  whereas  a  typical 
electron  current  for  this  pressure  range  in  N2  is  2  kA.  It  is  readily  apparent 
that  the  ion/electron  current  ratio  is  approximately  /(m^/M^^)  and  that  ion 
current  is  not  enhanced  appreciably  in  this  diode.  However,  if  the  charge 
exchange  or  recombination  processes  are  Included  in  the  above  estimate,  it  may 
be  found  that  the  ion  current  is  enhanced. 
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Many  Thomson  spectrograms,  similar  to  the  one  shovm  in  Figure  5-1,  reveal  a 
constant  peak  electric  deflection  angle  for  two  or  more  parabolas. This  is 
attributed  to  the  diode  like  accelerating  electric  field  which  results  in  a 
constant  peak  energy  per  charge  or  electric  deflection  as  is  readily  apparent 
from  Equation  (4-6).  Unfortunately,  this  fact  is  not  clear  in  most  spectrograms 
which  seem  to  show  a  constant  peak  velocity.  This  constant  velocity  line, 
however,  may  be  attributed  to  the  phenomena  of  charge  exchange  or  recombination 
which  may  occur  as  the  ions  traverse  a  high  density  plasma  and  a  neutral  gas. 

The  manner  in  which  this  occurs  is  illustrated  in  the  schematic  spectrogram 
shown  in  Figure  5-2.  All  of  the  ions  may  be  initially  accelerated  to  a  constant 
energy  per  charge  indicated  by  the  vertical  line.  As  the  charge  of  t''e  peak 
energy  ions  is  lowered  by  some  process,  they  will  be  recorded  on  the  next  lower 
parabola  with  the  charge  reduced  by  one  but  with  the  same  mass,  along  a  constant 
velocity  line.  In  this  way  ions  which  are  accelerated  while  they  are  in  a 
higher  charge  state  by  the  diode  accelerating  voltage  may  obtain  a  higher  energy 
per  charge  by  the  time  they  reach  the  spectrometer.  This  is  probably  the  source 
of  the  Mev  per  charge  ions  recorded  on  some  spectrograms.  Note  that  in 
Figure  4-8,  the  peak  proton  energy  is  about  200  keV  which  is  comparable  to  the 
peak  electron  energy  and  also  the  transmission  line  voltage  probe  peak  voltage. 

COMPARISONS 

Table  5-1  shows  results  of  experiments  done  at  other  laboratories  in  the 

past  15  years  on  the  charged  particle  beams  produced  by  a  plasma  focus.  Also 

21 

Figure  5-3  (due  to  Stygar  )  shows  the  scaling  law  for  electron  beam  current 
versus  bank  current  for  plasma  focus  devices  along  with  the  result  from  this 
experiment.  In  the  case  of  Reference  21,  the  diagnostic  was  similar  to  the  one 
in  this  work.  The  lack  of  agreement  is  perhaps  an  indication  of  the  fact  that 
the  experiments  performed  in  this  work  involve  a  different  filling  gas  pressure 
regime  which  may  be  more  favorable  to  beam  production. 
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TABLE 

5-1.  COMPARISON  OF 

MAJOR  CHARGED 

PARTICLE 

BEAM  RESULTS 

LAB 

DEVICE  PARTICLE 

DIAGNOSTIC 

ENERGY 

(keV) 

NUMBER 

Livermore 

900kA, 

2TorrD 

d 

nucl . act . 

>330 

ioi7 

800kA 

ITorrD 

e 

x-ray 

ave-150 

1.2kJ 

Darmstadt 

220kA 

4TorrD 

d 

magn.spect. 

80-400 

? 

220kA 

8TorrH 

e 

x-ray 

100 

? 

Stutt . 

1 

520kA 

-lOTD/Ar 

d,C,N,0  Thomson 

>350 

? 

Illinois 

560kA 

3TorrD 

d 

TOF/F.Cup 

25-lk 

10^7 

e 

Rogowski 
magn . defl . 

10-400 

17kA 

Stevens 

5-800kA 
3-5T  D 

d 

TOF/damage 

300 

10^^ 

e 

Rog . /dendrite 

s  300 

5kA 

NASA 

IMA 

STorrD 

e 

x-ray 

F.  cup 

<500 

lOkA 

This  Woi'k 

200kA 

.1-2T 

P.C.N, 
He ,  Ar 

Thomson 

20-2k 

10^^ 

e 

F.cup/Rog. 

magn.spect. 

20-300 

0.5-8  kA 

5- 
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CONCLUSIONS 

A  physical  model  for  the  plasma  focus  device  was  investigated.  The  diode 
model,  or  stationary  electric  field  model  of  the  plasma  focus,  has  been 
supported  to  the  greatest  extent  by  the  Thomson  spectrometer  diagnostic.  More 
support  for  the  diode  model  arises  from  the  fact  that  the  transmission  line 
voltage  probe  is  in  agreement  with  both  the  electron  spectrometer  and  Thomson 
spectrometer  peak  energies.  These  diagnostics,  however,  are  relatively 
insensitive  to  the  high  energy  tail  of  both  electron  and  ion  beams  which  is 
observed  at  many  laboratories.  The  yield  of  this  tail  is  several  orders  of 
magnitude  below  the  yield  in  the  range  investigated  here.  It  is  possible  that 
the  origin  of  these  particles  may  be  different. 

Several  new  diagnostics  were  created  for  this  experiment.  A  compact 
magnetic  electron  energy  analyzer  was  used  to  determine  the  energy  of  the 
electron  beam  in  the  plasma  focus  device.  The  resolution  and  trajectory 
analysis  for  this  diagnostic  was  performed  showing  the  resolution  to  be  better 
than  2  percent  at  the  lowest  energy  range.  The  observed  electronic  energies 
ranged  up  to  and  above  200  keV  in  a  device  with  charging  voltages  of  10  to  20 
kV,  a  better  than  10:1  stepup.  In  addition,  the  total  e-beam  energy  was 
measured  and  found  to  be  a  few  tens  of  Joules  with  an  efficiency  of 
approximately  1  percent.  The  quality  of  the  plasma  focus  produced  electron  beam 
was  investigated  for  the  first  time  in  this  work.  The  electron  beam  rms 
emittance  is  found  at  the  exit  of  the  center  electrode  drift  tube.  This  may  be 
the  place  where  the  beam  may  be  extracted  for  any  application. 

The  Thomson  spectrometer  was  used  to  obtain  the  energy  spectra  of  five 
different  filling  gas  ions.  In  addition,  the  Thomson  spectrometer  was  discussed 
and  the  resolution  analysis  performed.  It  is  hypothesized  that  charge  exchange 
and/or  recombination  is  responsible  for  the  presence  of  neutral  particles 
appearing  at  the  origin  of  the  Thomson  parabola,  as  well  as  the  presence  of  a 
constant  velocity  line  at  the  peak  energies  in  many  parabolas.  The  ion  beam 
flux  is  measured  with  a  CR-39  nuclear  track  detector  and  found  to  be  much  lower 
than  observed  by  other  researchers  in  devices  of  slightly  higher  currents.  An 
estimate  of  the  flux  was  also  made  with  the  results  of  the  ion  emittance  meter 
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and  was  found  to  be  consistent  with  the  result  of  the  fluxmeter.  The  ion  beam 
rms  emittance  is  investigated. 

This  work  is  a  comprehensive  evaluation  of  the  beam  behavior  of  the  plasma 

focus  produced  particle  beams.  In  the  future,  emphasis  should  be  placed  on  the 

use  of  time  resolved  diagnostics  which  may  yield  much  more  information  about  the 

nature  of  particle  acceleration.  In  particular,  time  resolved  energy 

measurements  of  the  electron  and  ion  beams  should  be  possible  through  the  use  of 

95 

new,  highly  sensitive,  streak  camera  systems.  In  conjunction  with  a  memory 
video  camera,  this  system  may  have  a  large  impact  on  data  acquisition. 

The  instabilities  that  occur  in  the  plasma  focus  subsequent  to  the  m-0  are 
not  yet  understood.  The  presence  of  microinstabilities  is  not  yet  well 
documented.  Future  investigations  should  be  placed  on  the  understanding  and 
controlling  of  the  rapid  rise  in  plasma  resistivity  which  leads  to  the  diode - 
like  nature  of  the  plasma  focus  observed  in  this  work.  Investigation  into 
insulator  properties  may  also  be  quite  important  since,  even  if  there  is  a  large 
induced  voltage  in  the  focus,  there  should  be  no  current  restrike  across  the 
insulator  to  short  this  voltage.  An  ideal  insulator  would  provide  a  uniform 
current  sheath  initiation  and  then  prevent  restrike  during  the  focus.  Future 
efforts  should  concentrate  on  making  the  device  more  reproducible;  an  ultra  high 
vacuum^^  may  improve  reproducibility  by  the  absence  of  impurities  in  the  plasma 
focus  chamber. 

Plasma  focus  research  in  the  past  15  years  has,  for  the  most  part,  been 
concerned  with  the  phenomenon  of  neutron  production  and  x-ray  production  in 
order  to  utilize  the  plasma  focus  as  a  pulsed  source  of  these  radiations.  In 
conjunction  with  these  experiments,  it  was  discovered  that  much  of  the  phenomena 
associated  with  the  plasma  focus  may  be  attributed  to  the  intense  charged 
particle  beams.  Much  progress  has  been  made  in  the  processes  involving  the 
beams;  however,  the  physical  mechanism  for  beam  production  is  still  a  mystery. 

It  is  hoped  that  this  research  will  aid  in  the  understanding  of  this  subject. 
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FIGURE  5-1.  A  THOMSON  SPECTROGRAM  WITH  A  CONSTANT  ELECTRIC 
DEFLECTION  SAME  VALUES  AS  FIGURE  4-7 
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APPENDIX  A 

MAGNETIC  ELECTRON  ENERGY  ANALYZER  RESOLUTION 

Resolution^' ^  is  an  important  consideration  in  the  design  of  an  energy 
analyzer.  For  convenience,  we  define  the  relative  energy  resolution  as  the 
uncertainty  in  the  measurement  at  some  energy  as  a  fraction  of  that 
energy,  AT/T  -  R.^,.  Under  the  assumption  of  uniform  magnetic  field  configuration 
provided  by  the  iron  case  and  magnet  as  well  as  the  absence  of  space  charge 
effects  in  the  analyzer,  the  resolution  is  determined  by  the  non- ideal 
collimation  of  the  two-slit  collimator.  Figure  A-1  shows  the  idealized  two- 
dimensional  geometry  required  for  resolution  analysis.  Three  different  particle 
trajectories  with  different  energies  all  impacting  at  the  same  point  on  the 
detector  are  shown.  The  detector  is  placed  along  the  line  (detector  line) 
separating  the  uniform  field  region  from  the  collimation  region.  In  this  case, 
the  downstream  collimator  slit  lies  directly  on  the  detector  line.  For  analysis 
we  define  new  coordinates  d  and  d,  where  d  is  the  distance  between  the  point 
where  any  arbitrary  particle  trajectory  crosses  the  detector  line  and  the  center 
of  the  downstream  slit  (see  Figure  A-1).  The  quantity  6  is  the  angle  that  the 
particle  trajectory  makes  with  a  line  normal  to  the  detector  line  shown  in  the 
figure.  If  the  limits  on  8  are  determined  by  a  double  slit  collimator,  then  the 
maximum  possible  8  is  just  half  the  maximum  acceptance  angle  of  the  collimator. 
If  the  upstream  and  downstream  slits  of  width  and  D2 ,  respectively,  are 

^Schneider ,  R.  F. ,  Luo,  C.  M.,  Rhee,  M.  J.,  and  Smith,  J.  R. ,  Rev.  Sci. 

Instrum.  56,  Aug.  1985. 
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separated  by  a  distance  L,  then  is  given  by: 


-1  ^1 

?  -  tan  "■{  - - -  )  . 

max  2L 


(A-1) 


Using  these  quantities  and  assuming  the  particles  follow  circular 
trajectories  upon  entering  the  B  field,  we  find  the  distance  from  the  center  of 
the  entrance  slit  to  the  point  of  impact  on  the  detector  as: 


z  “  2p(T)costf  +  d  , 


(A-2) 


where ; 


^  4-  9Tm,'^V^2 


CT  +  2Tmc 
ecB 


(A-3) 


The  equation  for  z  can  be  expanded  in  series  around  ^-0,  d  -  0,  and  T  -  T. 
Keeping  only  the  lowest  nonvanishing  order  terms,  we  arrive  at: 


2 

3z  1  d  z  dz 


(A-4) 


Suppose  we  measure  the  distance  z  exactly  (Az-0) .  The  energy  of  the 
electron  at  that  point  may  be  associated  with  that  corresponding  to  an  orbit 
with  radius  z/2 .  We  need  to  know  the  maximum  uncertainty  introduced  into  T 


(i.e.,  AT  )  by  A6  and  Ad  .  If  we  have  finite  Aff  and  Ad  ,  then  AT  will  iust 
max  max 

be  the  sum  of  the  magnitude  of  the  uncertainties  due  to  LS  and  Ad 

max  max 

Solving  for  AT  and  setting  Ad  -  D_  /2  and  Aff  -8  ,  the  resolution 

®  max  2  max  max 

function  is  found  to  be: 
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Figure  A- 2  is  a  plot  of  this  function  for  the  values  -  125/im  ,  D2  “ 
25/iin  ,  L  -  35  mm,  and  B  -  0.53  Tesla.  These  are  the  values  used  in  our 
experiment. 
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FIGURE  A-2.  RESOLUTION  FUNCTION  FOR  THE  MAGNETIC  ANALYZER 
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APPENDIX  B 

THOMSON  SPECTROMETER  RESOLUTION®*^ 

A  major  purpose  of  the  Thomson  spectrometer  is  to  separate  parabolas  of 
different  charge -to-mass  ratio  in  order  to  differentiate  ion  species  and  charge 
state.  This  in  principle  happens  naturally  if  the  parabola  is  an  infinitely 
thin  line.  Practically,  however,  this  can  never  be  done  since  the  collimation 
procedure  utilizes  two  pinholes  which  produce  an  image  with  finite  diameter  as 
displayed  in  Figure  B-1.  The  resulting  widths  of  the  parabolas,  which  are  swept 
out  by  the  different  energies  of  a  given  charge- to-mass  ratio,  leads  to 
overlapping,  which  reduces  the  resolution  of  the  spectrometer.  From  the 
geometry  of  Figure  B-1,  the  width  of  the  parabola  on  the  detector  is  given  by: 

^2 

"  L  ^°1  °2  • 

This  corresponds  to  a  subtended  angle  of: 

®3 

6  -  ■  (B-2) 

2 

The  overlapping  of  the  parabolas  as  they  near  the  origin  can  be  calculated, 
and  resolution  criterion  can  be  obtained  by  determining  the  points  where 


®' ^Schneider ,  R.  F. ,  Luo,  C.  M. ,  and  Rhee,  M.  J.,  J.  Appl .  Phys.  57, 
1985,  p.  1. 
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different  charge-to-mass  ratio  parabolas  are  separated  for  either  constant 
charge  or  constant  mass.  It  is  interesting  to  note  that  these  equations  arise 
solely  from  the  pinhole  geometry  and  are  independent  of  the  location  of  the 
field  region  so  that  they  are  applicable  to  the  conventional  Thomson 
spectrometer.  The  only  difference  between  compact  and  conventional 
spectrometers  is  the  location  of  the  pinholes  relative  to  the  fields.  In  the 
compact  spectrometer,  the  field  region  is  placed  between  the  first  and  second 
pinholes  allowing  much  more  control  over  the  pinhole  collimation  image  on  the 
detector.  In  conventional  spectrometers,  the  pinhole  collimation  system  is 
entirely  upstream  of  the  fields.  This  implies  the  distance  L2  in  Figure  B-1 
must  be  large  enough  to  contain  the  fields.  This  limits  the  minimiim  obtainable 
and,  hence,  requires  a  spectrometer  which  is  an  order  of  magnitude  larger  to 
obtain  the  same  resolution.  There  is  a  slight  condition  which  must  be  met  in 
the  use  of  this  spectrometer,  namely  the  entire  spectrometer  entrance  must  be 
uniformly  illuminated  by  the  ion  source.  This  is  quite  easy  to  achieve  in  the 
plasma  focus  since  the  source  is  of  finite  spatial  dimension  and  the  first 
pinhole  is  40  cm  upstream  of  the  spectrometer  and  8  cm  upstream  of  the  focus 
allowing  a  source  of  less  than  0.1  mm  to  cover  the  entire  entrance  to  the 
spectrometer . 

To  analyze  charge  resolution,  we  suppose  that  we  have  one  ion  species  with 
several  different  charge  states.  Along  what  curve  can  we  find  points  at  which 
neighboring  parabolas  are  separated?  To  answer,  we  use  the  criterion  for 
separation  that  the  separation  of  the  centers  of  the  parabolas  along  a 

vertical  line,  i.e.,  6^  direction  must  be  greater  than  the  parabola  width  along 
that  line: 
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hB  >  6(1 
m 


6  1/ 

+  -a. 

hi 

e 


and  LB 


(B-3) 


This  is  illustrated  in  Figure  B-2.  Here  we  have  made  use  of  the  fact  that 
the  slope  of  the  parabola  is  given  by; 


dB  B 
_ m  _ m 

dB  “  IB  ’ 
e  e 


(B-4) 


which  is  derived  from  the  parabola  equation.  Assuming  small  we  can 

substitute  the  separation  condition  above  into  the  total  differential  of  the 
parabola  equation  with  Z,  A,  6^,  and  6^  as  variables  yielding: 


u  lEdl  ^^m  /JI!_n2,V2 

'  (jBdl,2  «e  ‘2.  >  1  • 


If  dA  -  0  and  dZ  -  1 ,  we  will  arrive  at  the  equation  describing  the  line 
connecting  the  points  of  charge  resolution  for  a  certain  ion  species: 


4  2  2 

m  m  e 


-  0 


(B-6) 


We  have  made  the  definition: 


k 


e  (iBdlJ^ 
u  JEdl 


The  physical  solution  to  the  quartic  equation  is: 


B 

e 


(B-7) 


(B-8) 


B-3 


NSWC  TR  88-394 


which  describes  a  straight  line  which  intersects  the  origin.  This  is 
illustrated  in  Figure  B-3.  Below  this  line,  the  parabolas  of  every  charge  state 
of  ion  mass  A  are  separated.  An  alternate  way  of  viewing  this  is  to  note  that 
such  lines  are  constant  velocity  lines  with  velocity: 

Any  ions  with  velocity  less  than  this  value  will  be  resolved.  Thus,  large 
values  of  k/6  will  yield  better  charge  resolution. 

For  mass  resolution,  let  dZ  -  0,  dA  -  1.  Using  the  parabola  Equation  (4-8) 
to  eliminate  A,  we  obtain: 

8^-  (ZkSe  -  4(Zk5fl  )^  -  0  ,  (B-10) 

mm  e 

which  can  be  solved  numerically.  An  interesting  analytical  result  can  be 
derived  when  the  charge  state  Z  in  the  equation  for  mass  resolution  is 
eliminated  in  favor  of  mass  number  A  through  the  parabola  equation.  In  this 
case,  we  find  the  curve  which  connects  the  resolution  points  (points  of 
separation)  of  parabolas  with  atomic  numbers  A  and  A+1  with  different  charge 
states : 


8 

m 


2ASe 


(B-11) 


This  equation  describes  a  hyperbola  with  asymptotes  8 ^-A6  and  6 


2AS  above 


which  isotopes  are  resolved  as  in  Figure  B-4.  This  can  be  generalized  to  dA  -  n 
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provided  n/A  is  much  less  than  1  in  order  to  obtain  a  resolution  curve  for 
elements  which  differ  in  mass  number  by  n.  This  leads  to  the  simple 
modification  in  the  previous  equation  that  A  must  be  replaced  by  A/n  everywhere 
it  appears. 


Energy  and  momentum  resolution  is  also  limited  by  the  subtended  angle  of 
the  pinhole  collimation  spot,  6.  The  spectrogram  can  be  thought  of  as  a 
superposition  of  many  of  these  spots  infinitely  close  together  along  the 
parabola  length.  Energy  or  momentum  which  is  found  from  either  the  electric  or 
magnetic  deflection  angle  is,  hence,  uncertain  by  the  an  amount  corresponding  to 
the  projection  of  the  collimation  spot  length  delta  along  the  parabola  onto  the 
electric  or  magnetic  deflection  axis  of  the  spectrogram  (see  Figure  B-5). 

Within  our  approximation,  these  quantities  are  related  by  the  equation; 

fi  -  [l+(^)  .  (B-12) 

e 


where  we  have  used  Equation  (B-4) .  Relative  energy  and  momentum  uncertainties 

AT/T  and  Ap/p  can  be  found  from  the  deflection  equations  in  terms  of  S  as , 

e 


a 


T 


T 


(B-13) 


(B-14) 


A  quantity  representing  either  T  or  p  can  be  expressed  along  with  its 
uncertainty  as  Q^(l  ~  aj/2).  It  is  of  interest  to  find  the  curves  of  constant 
which  will  give  a  contour  map  of  relative  uncertainty  of  energy  or  momentum. 
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Combining  the  equations  for  6  and  a  yields: 


- - -  +  - - - 1 

(5/0^)^  (26/q^)^ 


_  e  m  , 

2  +  2  ”  ^  ■ 

(5/2Qp)^  <VQp) 


(B-15) 


(B-16) 


These  describe  upright  ellipse  about  the  origin  which  are  contours  of 
constant  relative  uncertainty  as  illustrated  in  Figure  B-6.  The  relative  energy 
or  momentum  resolution  at  a  certain  point  on  the  parabola  corresponds  to 
the  or  Qp  value  of  the  ellipse  that  intersects  that  point.  In  terms  of 
spectrometer  parameters,  we  may  find  a.^.  by  combining  the  deflection  Equations 
(4-6)  and  (4-7)  with  Equation  (B-16). 


S_  ^  Z  .  _k_  ,  ,  -  V2 

®t‘  «  A 

e  e 


(B-17) 


The  resolution  function  is  plotted  in  Figure  B-7  for  the  spectrometer  parameters 
used  for  the  energy  spectrum  measurements. 
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FIGURE  B-2.  SEPARATION  OF  PARABOLAS  FOR  CHARGE  RESOLUTION  ANALYSIS 
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FIGURE  B-3.  CHARGE  STATE  RESOLUTION  OF  THE  THOMSON  SPECTROMETER. 

THE  LINE  SEPARATES  THE  REGIONS  WHERE  CHARGE  IS  RESOLVED 
PROM  REGIONS  WHERE  IT  IS  UNRESOLVED. 
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FIGURE  B-6.  CONTOURS  OF  CONSTANT  RELATIVE  ENERGY 
OR  MOMENTUM  UNCERTAINTY  WHICH  ARE 
ELLIPSES 
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FIGURE  B- 


ENERGY  PER  CHARGE  T/Z(MeV) 


PLOT  OF  RELATIVE  UNCERTAINTY  VERSUS  ENERGY  PER  CHARGE 
FOR  THE  SPECTROMETER  USED  TO  OBTAIN  ION  ENERGY  SPECTRA 
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